










« v'' v ^l< ' •'' '.vv \ "* $ !sM«$5 m m is 'tm s§^^ss3& &§£&& 


' -r': 

■N0 


Ki 






; 


National Ae. 
||paca Adml 


X95— 10160 


( NASA— CR— 4636 ) EXPERIMENTAL WATER 
DROPLET IMPINGEMENT DATA ON 
AIRFOILS, SIMULATED ICE SHAPES , AN 
ENGINE INLET AND A FINITE WING 
Final Report (Wichita State Univ.) 
178 p LIMIT USGA AND C 


Unci as 


El/03 0035344 




NASA Contractor Report 4636 
DOT/FAA/CT-TN93/1 8 


Experimental Water Droplet Impingement Data on 
Airfoils, Simulated Ice Shapes, an Engine Inlet 
and a Finite Wing 


M. Papadakis and X. Liu 
Wichita State University 
Wichita, Kansas 

M. Breer and N. Craig 
Boeing Military Airplanes 
Wichita, Kansas 


Prepared for 

Lewis Research Center 

under Contract NAG3-566 



National Aeronautics and 
Space Administration 

Office of Management 

Scientific and Technical 
Information Program 

1994 



U.S Department 
of Transportation 

Federal Aviation 
Administration 


DOCUMENT AVAILABILITY CHANGE NOTICE 


NASA Contractor Report 4636 


Experimental Water Droplet Impingement Data on Airfoils, 
Simulated Ice Shapes, an Engine Inlet and a Finite Wing 

M. Papadakis, M. Breer, N. Craig, and X. Liu 


December 1994 


This document was published in December 1 994 
with a restriction. 

It was changed March 2005 to Unclassified/Unlimited per DAA modified March 30, 2005. 


Per the STI Program Office and Code I at HQ, 
you may modify copies in your possession. The restriction notice 
on the cover, title page and RDP should be boldly crossed out 
and the above statement printed clearly above or below it. 


Glenn Research Center 
Publishing Services Coordination Office 
216-433-3207 


Experimental Water Droplet Impingement Data on Airfoils, 
Simulated Ice Shapes, an Engine Inlet and a Finite Wing 


M. Papadakis, M. Breer, N. Craig and X. Liu 


ABSTRACT 

An experimental method has been developed to determine the water droplet impinge- 
ment characteristics on two- and three-dimensional aircraft surfaces. The experimental 
water droplet impingement data are used to validate particle trajectory analysis codes 
that are used in aircraft icing analyses and engine inlet particle separator analyses. The 
aircraft surface is covered with thin strips of blotter paper in areas of interest. The surface 
is then exposed to an airstream that contains a dyed- water spray cloud. The water droplet 
impingement data are extracted from the dyed blotter paper strips by measuring the opti- 
cal reflectance of each strip with an automated reflectometer. Experimental impingement 
efficiency data are presented for a NLF(l)-0414 airfoil, a swept MS(l)-0317 airfoil, a Boe- 
ing 737-300 engine inlet model, two simulated ice shapes and a swept NACA 0012 wingtip. 
Analytical impingement efficiency data are also presented for the NLF(l)-0414 airfoil and 
the Boeing 737-300 engine inlet model. 
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EXECUTIVE SUMMARY 


The water droplet impingement characteristics of an aircraft are typically used to 
determine size and location of ice protection systems. The impingement characteristics 
have been historically determined experimentally. With the advent of powerful computers, 
a number of computers codes have been developed capable of predicting the impingement 
characteristics of aircraft surfaces and performing icing analyses. To insure that these codes 
are correct, they must be validated against experimental data that must span the range 
of geometries and flight conditions to which the computer codes will be applied. Until 
recently, the experimental data base was limited to two-dimensional cylinders, airfoils, 
several bodies of revolution and a supersonic inlet. These data were generated by the 
National Advisory Committee for Aeronautics (NACA) in the 1950’s using a dye tracer 
technique and colorimetric analysis. 

To provide a more extensive water droplet impingement database, NASA and the FAA 
began a research program in August 1984 with Wichita State University and Boeing Mil- 
itary Airplanes. The main objectives and accomplishments of this work were presented in 
NASA CR 4257 published in November 1989. Briefly, experimental and data reduction 
methods were developed and verified for measuring local droplet impingement efficiency 
on a variety of aerodynamic surfaces. The data reduction method was based on laser re- 
flectance which is more accurate and considerably more efficient than colorimetric analysis. 
An automated data reduction system was also developed to make practical the reduction 
of a large mass of data. Water droplet impingement data were obtained for a 4-inch di- 
ameter cylinder, a NACA 652015 airfoil, a supercritical airfoil, an axisyxnmetric inlet and 
a Boeing 737-300 engine inlet. Experimental data were correlated with analytical results. 

The work presented in this report is the conclusion of the research program which was 
initiated in 1984. The following tasks were accomplished: 

1. The design of the collector mechanism for measuring LWC was improved. 

2. Water droplet impingement data were obtained for a Natural Laminar Flow 
airfoil, a swept MS-0317 airfoil, a swept wingtip, a Boeing 737-300 engine 
inlet and two simulated ice shapes. 

3. Experimental data and computational impingement results were correlated 
for selected geometries and test conditions. 
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1.0 INTRODUCTION 


Aircraft flying at subsonic speeds below about 8000 meters are subject to ice accretion 
caused by very small supercooled water droplets that freeze upon impact with the aircraft. 
Ice accretion, which occurs on the forward facing surfaces of an aircraft, can be detrimental 
when occuring on the leading edge of a wing and/or on engine inlet diffusers. On wing 
surfaces, ice accretion can lead to an increase in drag, resulting in excessive fuel consump- 
tion or flow separation, which results in a decrease in lift. Engine inlet diffusers are more 
susceptible to the effects of ice accretion, particularly when shedding ice fragments enter 
the engine inlet and cause mechanical damage to the turbomachinery of the engine. Thus, 
it is necessary that ice protection systems, for the removal or prevention of ice accretion, 
are incorporated in aircraft. 

Current ice protection systems can be divided into two different categories: de-icing sys- 
tems and anti-icing systems. De-icing systems periodically shed ice that has already formed 
by mechanical, thermal or chemical means. Deicing systems have low energy requirements 
but have higher aerodynamic penalties. Anti-icing systems prevent ice formation by means 
of thermal energy or a Freezing Point Depressant (FPD). While aerodynamic penalties are 
held to a minimum with an anti-icing system, the energy requirements are higher than 
those of a de-icing system. 

Of major concern in the design of an anti-icing or de-icing system is the extent and 
amount of water impingement. The Impingement characteristics of an aircraft can be used 
to determine size and location of ice protection systems. The impingement characteristics 
have been historically determined experimentally but recent advances in computer speed 
have allowed the use of trajectory codes to determine them analytically. 

With the advent of powerful computers, a number of computer codes have been devel- 
oped (Ref. 1-5) capable of predicting the impingement characteristics of aircraft surfaces 
and performing icing analyses. These codes can provide a cost effective means for the 
design of ice protection systems. To insure that these codes are correct, they must be 
validated against experiment. The experimental data must span the range of geometries 
and flight conditions to which the computer codes will be applied. Until recently, the data 
base was limited to two-dimensional cylinders, airfoils, several bodies of revolution and a 
supersonic inlet. These geometries were tested by NACA in the 1950’s using a dye-tracer 
technique (Ref. 6-9). 

NASA and the FA A began a research program in August of 1984 with Wichita State 
University and Boeing Military Airplanes to provide a more extensive water droplet im- 
pingement data base. The objectives of this program are listed below: 


a . Develop an efficient and accurate method for measuring local droplet im- 
pingement efficiency on a variety of aerodynamic surfaces in two and three 
dimensional flows. 

b . Develop experiments to verify the method. 
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c . Test several state-of-the-art airfoil sections as well as wing and engine inlet 
geometries to produce a comprehensive experimental water droplet impinge- 
ment data base. 

d . Develop an automated data reduction system that would be practical and 
efficient for the acquisition and reduction of large quantities of data. 

Several two and three dimensional models were tested during the initial phases of this 
research program (Ref. 10) . The current tests represent the completion of data acquisition 
for the water droplet impingement data base. This report presents the testing of several 
modern geometries: a Natural Laminar Flow airfoil, a swept MS-0317 airfoil, a swept 
wingtip, a Boeing 737-300 engine inlet and two simulated ice shapes. 
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2.0 TRAJECTORY EQUATION AND IMPINGEMENT PARAMETERS 

The non-dimensional droplet trajectory equation and impingement parameters which 
are commonly used in the presentation of theoretical and experimental impingement data 
are presented in this section, A more detailed discussion of impingement parameters can 
be found in Ref, 10, 

2.1 Equations of Particle Motion 

The forces on a small spherical droplet moving in the steady flow of air are droplet 
drag, weight and buoyancy. The predominant force acting on a droplet is the fluid dynamic 
drag arising from the relative (slip) velocity of air with respect to the droplet (Ref. 10). 
The particle equation of motion in non-dimensional form is: 

dUj _ - Vi) (1 - a) gLSq 

dt 24 K ^ 

where 

— Cd = drag coefficient of spherical particle 

— K — PparUcicUood 2 inertia parameter of droplet 

— t = time, dimensionless with L/Voo 

— <7 = p/ppy density ratio of air to particle 

— L = characteristic dimension of body 

— R v — Pair\V — relative Reynolds number of droplet 

— Ui — *th component of particle velocity, non-dimensionalized by 

— Vi = tth component of _air velocity, non-dimensionalized by 

— &ij = Kronecker 6 

— d — particle diameter 

— p = absolute air viscosity 

— p = air density 

— p p = particle density 

— Voo = freestream air velocity 

* = 1 x — direction , i = 2 y — direction , * =3 z — direction 
The derivation of equation (2 — 1) is based on the following assumptions: 

• Single phase (air) flow about the body (i.e., air flow is not disturbed by the 
presence of the droplets) 
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• Quasi-steady-state approximation (i.e., at each instant and position the 
steady-state drag and other forces act on the particle) 

• Compressible or incompressible potential flow field of the gas phase about 
the body 

• Particles of spherical shape 

• Viscous flow effects such as thick boundary layer formation and flow separa- 
tion are not considered 

The above mathematical model is a valid approximation for typical icing conditions 
within the intermittent and continuous maximum icing envelopes of the PA A (Ref. 11). 
The maximum concentration and mean volumetric diameter (MVD) for these conditions 
are: 

Intermittent Maximum Continuous Maximum 

LWC Max 3.0 g/m B LWC Max 0.8 g/m B 

MVD Max 50 MVD Max 40 #tm 

For the concentrations and sizes of droplets expected to occur within icing clouds, the 
assumptions of undisturbed airflow and spherical shape (due to surface tension) of droplets 
are very closely approximated. 

Given an expression for Cd(R v ), the drag coefficient of the droplet as a function of the 
relative Reynolds number of the droplet, equation (2 — 1) can be integrated numerically 
to find the particle trajectory path. An example of Cd as a function of R v can be found 
in Ref. 10. 

2.2 Impingement Parameters 

Impingement parameters for non-uniform droplet distributions are discussed below. 

2.2.1 Liquid Water Content (LWC) 

The liquid water content (LWC) of a cloud is the amount of water contained in a given 
volume of cloud. LWC is usually expressed in grams of water per cubic meter of cloud. 
Typical LWC maI values for icing clouds are listed in section 2.1. 

For simulated icing clouds, the LWC is controlled by the water and /or air pressures 
of the spray system used to create the spray clouds. 

2.2.2 Cloud Droplet Distribution 

The distribution of droplets in a cloud can be expressed in many various forms (Ref. 
10), but for convenience and computational efficiency, a discrete form of cumulative liquid 
water content versus droplet diameter normalized to mean volumetric diameter is used. 
The droplet distributions used in this work were obtained experimentally as explained in 
Ref. 10 and are given in Table 2.1 and 2.2. 

A distribution which has been employed in various analytical studies is the Langmuir 
“D” . This distribution and other similar distributions were established by Langmuir by 
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matching the data from natural-icing cloud measurements made on Mt. Washington (Eef. 
12). Comparisons of Langmuir “D* distribution and droplet distributions produced by the 
spray system used in the current tests are shown in Fig. 2.1. 

2.2.3 Mean Volumetric Diameter (MVD) 

By definition, the Mean Volumetric Diameter (MVD) of a droplet distribution is that 
mean diameter for which half the total liquid water content is contained in droplets larger 
than the mean and half in droplets smaller than the mean. Given a droplet distribution 
the MVD van be calculated as follows. 

For a continuous distribution if n(D) is the number of particles per unit sampling 
volume having diameters between D and D+dD (volumes between V and V+dV) then 
Dmv d can be calculated from 


If 

5 e* IS::; n(z)x'dx 


= 0.5 


(2 - 2 ) 


For a discrete distribution if the particle number density is given in N discrete groups 
such that n t (Di) is the number of the particles in group i having diameters between D and 
D+dD then, equation (2 - 2) can be written as 


s*. e mol 

=0.5 (2-3) 

I p . E 

*=1 


where: 

Djc = T*mvd 

p w ~ density of water, g/cm 3 


2.2.4 Local Impingement Efficiency (0 ) 

Considering a body in a cloud with uniform droplet size distribution the local impinge- 
ment efficiency 0 for any point on the body surface is defined as the local droplet flux rate 
at the body surface normalized to the free stream flux rate. Referring to figure 2.2 0 is the 
ratio of that infinitesimal area dAoo to the corresponding impingement area on the body 
surface dA a . This definition follows from the continuity of droplet mass flow. 

For a continuous non-uniform cloud distribution, the impingement efficiency is given 
by equation (2-4). 

0=1-1 0dw (2-4) 
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where 0 is a function of drop size and therefore can be expressed as a function of w, the 
liquid content for a given drop size. 

For a discrete cloud distribution such as the one shown in Fig. 2.1, 0 is defined as the 
weighted average of the local impingement efficiency values due to each droplet group in 
the cloud. Let w t be the liquid water content of the cloud, A Wi be the partial liquid water 
content contained in the droplets of size (d»), in the group (i) of the distribution, and N be 
the total number of discrete size droplet groups available. For a body exposed to a cloud 
with such a droplet distribution, the local impingement efficiency due to a single droplet 
group of size d, is where 0 is defined in Fig. 2.2. The local impingement efficiency 
due to all N groups in the distribution over an infinitesimal area of the body is given by 
equation (2 — 5) below. 

1 N 

(2-5) 

* t-i 

2.2.5 Total Impingement Efficiency ( E ) 

The total impingement efficiency of a three dimensional body exposed to a cloud of 
droplet distribution is 



where: 

A/ is the projected frontal area of the body 

dA t is an infinitesimal impingement area on the surface of the body 

In order to integrate equation (2—6), 0 must be known as a function of surface location. 
Such a function can be defined from experimental or analytical results. 

2.2.6 Impingement Limits ($ m /L) max 

Droplets which start out at a freestream position (t/oo) with respect to a reference line 
that pass through the highlight ( most forward point at a = 0° ) of a body downstream 
will impinge at some position on that body. As these initial freestream droplet positions 
increase in distance from the reference line they will impinge farther back along the surface 
of the body until a maximum distance (j/oo, max) is obtained. This limiting trajectory is 
defined as the tangent trajectory to the body at point P (Fig. 2.2). Any droplet starting 
at a freestream position farther from the reference line than (t/oo, max) will miss the body 
entirely. The distance s m measured along the body surface from the highlight of the body 
to the point P is called the limit of impingement. This distance is usually expressed in 
dimensionless form by dividing s m by the characteristic length (L) of the body. 

For two-dimensional flow, there are two impingement limits, an upper and a lower (for 
external flow, e.g., airfoil section) or an outer and inner ( for partly internal flow, e.g., 
engine inlet). For three-dimensional flow, the limits of impingement may vary spanwise 
along the surface of a finite wing or circumferentially along the surface of an engine inlet. 
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For a droplet distribution that varies from d min to d mazi the impingement limits can 
be established for each droplet size. The maximum impingement limits are defined by the 
impingement limits of the largest droplet diameter of the distribution. 

2.2.7 Local Impingement Rate 

To calculate the amount of water collected on a small area of the body where the local 
impingement efficiency value is known, the following equation can be used. 

Wp = 0.3296 Voottf*/? (2-7) 

where: 

Wp = Local impingement rate, lb/ hr ft 2 
Voo = frees tream velocity, mph 
wt = Total liquid water content, gjm 3 
ft = Local impingement efficiency 
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TABLE 2.1 

DISCRETE CUMULATIVE DISTRIBUTION 
FOR 0.80 PRESSURE RATIO 


CUMULATIVE % LWC(I) 

% LWC(I) 

D(I)/MVD 

97.5 

5 

0.3161 

90.0 

10 

0.4981 

75.0 

20 

0.6872 

50.0 ’ 

r ~” 30 

1.0000 

25.0 

20 

1.3737 

10.0 

10 

1.9614 

2.5 1 

5 

2.8288 


MVD = 16 MICRONS 


TABLE 2.2 

DISCRETE CUMULATIVE DISTRIBUTION 
FOR 0.65 PRESSURE RATIO 


CUMULATIVE % LWC(I) 

% LWC(I) 

D(I)/MVD 

97.5 

5 

0.2770 

90.0 

10 

0.4460 

75.0 

20 

0.6617 

50.0 

30 

1.0000 

25.0 

20 

1.5865 

* 10.0 

10 

2.2943 

2.5 

5 

3.2542 


MVD = 20 MICRONS 
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CUMULATIVE % LWC CUMULATIVE % LWC 



(A) MVD = 16 fxm 



(B) MVD = 20 pm 
FIGURE 2.1 

EXPERIMENTAL AND LANGMUIR D CUMULATIVE DROPLET DISTRIBUTIONS 
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FIGURE 2,3 

TWO-DIMENSIONAL DROPLET TRAJECTORIES FOR A BODY 
IN A CLOUD OF UNIFORM DROPLET SIZE 


3.0 EXPERIMENTAL METHOD 


3.1 Test Approach 

The dye-tracer technique that was employed in the water droplet impingement tests 
conducted in September of 1985 was used for the current tests. Details of the experimental 
approach are given in Ref. 10. Briefly, distilled water containing a known concentration 
of blue dye was injected into the airstream of a wind tunnel in the form of a droplet spray 
cloud. The test model was covered with thin strips of blotter paper in areas of interest 
and was exposed to the spray cloud. The amount of dye-mass per unit area of blotter 
strip obtained in a given time interval was measured using a laser reflectance method. The 
water impingement characteristics of a surface was then extracted from the concentration 
and location of the dye distribution on the blotter paper. 

To calculate the local impingement efficiency from the measured dye-mass per unit 
area, an accurate measurement of the LWC in the spray cloud is required. In reality, a 
totally uniform spray cloud is very difficult to achieve. During the current tests LWC 
was measured at each location in the wind tunnel test section where a blotter strip was 
placed on a model. This was accomplished by using a reference collector mechanism with 
nine small collector blades. The impingement efficiency of these blades, for the droplet 
distributions given in tables 2.1 and 2.2, were close to 100 percent. The collector mechanism 
was placed in the test section of the wind tunnel with its collector blades positioned as 
close as possible to the blotter strip locations on the test model. A blotter strip was placed 
on each collector blade and impingement tests were conducted at identical conditions that 
the model was tested. During all collector tests the test model was removed from the test 
section. Using the laser reflectance method on the reference collector blotter strips the 
amount of dye per unit area of freestream at all locations of interest was determined. The 
local impingement efficiency, /? , at a given model location was obtained from equation 

( 3 - 1 ). 

Measured amount of dye/unit surface area/unit time on model 
at a given location 

0 = - (3 - X) 

Measured amount of dye/unit surface area/unit time in the 
freestream for a collector whose 0 = 1 

The reference collector that was used for the current tests is discussed in Section 3.5. 

It should be noted that the same blotter paper from the same batch was used for the 
April 1989 tests as it was in the September 1985 tests. This is also true for the blue dye 
used in these tests. This is important when validating the laser calibration curve that is 
used in the data reduction method. 
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3.2 Test Facility 

All of the water droplet impingement tests described in this report were conducted 
in the NASA Lewis Icing Research Tunnel (IRT). This tunnel has a test section of 6 ft 
high by 9ft wide and is capable of an airspeed of 300 mph with air temperatures as low 
as -30 °F. The IRT has a spray system that produces simulated icing clouds with LWC 
ranging from 0.5 to 2.0 g/m z . 

All impingement tests were conducted at a temperature of 35 to 70 °F (i.e., no icing) 
and at a true air speed of 165 mph. The IRT spray system was not used to produce the 
spray clouds for the impingement tests for reasons explained in Section 3.3. However, the 
IRT spray system was activated occasionally between test runs to help maintain a high 
humidity level in the tunnel airstream and keep droplet evaporation low. The spray system 
designed and used during the 1985 tests was employed in all 1989 impingement tests. This 
system is briefly described in section 3.3, details of the spray system design can be found 
in Ref. 10. 

3.3 Spray System 

To avoid saturation of the blotter strips with the dyed water, a short spray interval 
was required, typically 2.6 to 4.5 seconds, depending on cloud MVD size. For the larger 
MVD size shorter spray durations were used. The IRT spray system was not designed 
for such short spray intervals, so an alternate system was designed that could accurately 
control the spray time interval. 

The spray system used for the 1989 tests used the same principles and components that 
were used in the 1985 tests. The spray system consisted of a 9 gallon aluminum tank (which 
contained the dyed water solution), a pressure system and twelve nozzles (Fig. 3.1), Each 
nozzle assembly consisted of a solenoid valve that controlled the water supply, a nozzle 
housing that contained the atomizing air inlet, a pressure gauge or pressure transducer to 
monitor water pressure, and a NASA standard icing spray nozzle (Fig. 3.2). 

The dye solution, which was pressurized in the tank to approximately 106 psig, passed 
through rubber hoses to each individual nozzle assembly. The air supply to the nozzle 
assemblies was controlled by a single solenoid valve at the main air supply. During testing, 
the air supply was turned on a few seconds before the water supply to allow the air pressure 
to stabilize. Next the dyed water supply was turned on and a spray cloud was produced. 
The mean volumetric diameter (MVD) of the spray cloud was controlled by varying the 
air-to- water pressure ratio. Two air pressures were used in the current impingement tests. 
These were 65 psig and 80 psig corresponding to clouds with MVDs of 20 and 16 microns 
respectively. 

The spray system was mounted on the existing spray bars of the IRT spray system 
located in the plenum chamber of the tunnel. The nozzles were initially set up in the wind 
tunnel at the same position as they were in September 1985 in an effort to duplicate the 
spray cloud uniformity in those tests. After checking the uniformity of the spray cloud, 
the nozzles were relocated until a satisfactory degree of cloud uniformity was obtained (see 
Section 3.4). 
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Once uniformity was achieved, more extensive measurements of the nozzle locations 
were taken. Originally, only the nozzle locations with respect to the spray bars in the 
IRT were recorded. During uniformity testing it was determined that the spray cloud is 
sensitive to nozzle orientation with respect to the tunnel floor and tunnel walls. To obtain 
a better definition of nozzle placement in the IRT, additional measurements were taken. 
With the tunnel at zero airspeed and with the spray system air supply off, the nozzles 
were turned on to allow a pressurized stream of dyed water to flow from each nozzle. The 
stream of water (as opposed to a spray when the air was on) exited from the nozzles and 
impinged on the tunnel floor near the tunnel bellmouth. The coordinates of the points of 
impingement were recorded for each nozzle. These additional measurements will allow for 
a more accurate recreation of the spray cloud should any further testing be required. The 
nozzle locations with respect to the IRT spray bars and the tunnel floor are shown in Fig. 
3.3. 

3.4 Uniformity Test 

Cloud uniformity was measured using a 6 ft x 9 ft brass grid with horizontal and 
vertical increments 6 inches apart. Blotter strips were attached to this grid to cover an 
area of 2 ft x 3 ft centered about the tunnel longitudinal axis, as shown in Fig. 3.4 The 
tunnel was brought up to test speed and the blotters were sprayed. The dye distribution 
on each blotter ( 35 blotter strips were used during each test) was determined using the 
data reduction method, and the nozzles were adjusted to remove dark or light regions of 
dye distribution. 

During the 1985 uniformity test runs, all spray nozzles were activated during each uni- 
formity test and problem areas in the spray cloud were identified by measuring variations 
in the color of the blotter squares and their perspective locations. In this method it was 
difficult to determine which nozzles should be moved to correct any light or dark regions 
in the array of blotter squares. In the current uniformity tests each nozzle of the spray 
system was tested individually to determine its contribution to the spray cloud. This gave 
the opportunity to make a more educated guess on how the positions of nozzles should be 
adjusted to obtain a uniform spray cloud. 

It was also decided that a single 2 to 3 ft long 1.25 inch wide strip of blotter paper 
could be used to replace a row or column of blotter squares (used in the 1985 tests), 
allowing for a continuous measurement of dye deposit along the row or column instead of 
measurements at six inch intervals (see Fig, 3.5). To prevent excess flow blockage in the 
area of the uniformity grid, the test was conducted separately for the five rows of strips 
and the seven columns of strips. 

After a satisfactory cloud distribution was obtained (by visual check), the dye deposits 
on the blotter strips were measured using the laser reflectance method to verify that the 
deposits did not vary more than ±10% from the average. Further adjustments were then 
made to the spray system based on the results obtained from the laser reflectometer to 
meet *he ±10% uniformity requirement, A total of 59 uniformity runs were performed (see 
Appendix B). 
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3.5 Reference Collector Mechanism 

To calculate the local impingement efficiency, the local LWC of the free stream spray 
cloud had to be determined by use of a “ reference collector mechanism”. The reference 
collector mechanism consisted of ah airfoil shaped vertical beam with eight arms spaced 
at 45 degree intervals that extended in front of the beam (Fig. 3.6). A 4 inch long wedge 
shaped collector blade was mounted on the end of six of the arms, and a 18 inch long 
collector blade was mounted between the top and bottom arms of the reference collector 
mechanism. The collector blade was designed to hold a 0.2 inch wide x 1.5 to 2 inch 
long rectangular blotter strip in place while minimizing the disturbance to the flow in 
the proximity of the collector blade. This reference collector was designed to obtain the 
local LWC for the airfoils and the engine inlet tested, replacing the two separate reference 
collectors that were used in the 1985 tests. 

The reference collector was positioned so that the blotter strips on the collector blades 
were at the same position with respect to the test section as the blotter strips that were 
mounted on the test models. The reference collector was then exposed to a spray cloud 
under the same test conditions as the models. This test was repeated for all test conditions 
and all angles of attack tested. 

Corrections to the measured LWC were made based on computer predicted blade col- 
lection efficiencies. Even though the reference collector was designed to maximize collection 
efficiency, the impingement efficiency of the collector blades was not 100%. If the blade 
collection efficiency were 100% then the LWC at the blade would be a direct measure of 
the LWC in the tunnel. The predicted collector efficiency values for the spray clouds with 
MVD’s of 16 and 20 microns were 86% and 89% respectively. The LWC can be defined as 
the amount of dye per unit time of collector exposure divided by the theoretical collector 
impingement efficiency and the concentration of the dye solution: 

Mass of dye per unit time on collector blotter strip//? collector 
for a given MVD 

LWC = (3 - 2) 

Concentration of dye solution 


3.6 Mass Flow Measurements 

The flow field about an engine inlet is largely dependent on the mass flow through the 
inlet. To obtain a broad water droplet impingement data base, it was necessary to test 
the engine inlet at different mass flow rates as well as different angles of attack. 

To provide the variations in mass flow, the suction system of the IRT was utilized. 
The system consisted of a large vacuum chamber and associated vacuum pumps outside 
the tunnel that were connected to the tunnel test section by a 24 inch diameter pipe. The 
suction was regulated by a valve that could be opened or closed by an electrical switch that 
was located inside the control room of the IRT. The tunnel test section vacuum system 
was connected to an elbow arrangement that extended through the ceiling of the IRT test 
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section. The elbow was connected to the engine inlet model mount by means of a flexible 
rubber hose, which allowed for variations in angle of attack. 

The mass flow through the system was measured with a KURZ velocity and tempera- 
ture probe that was inserted in the 24 inch diameter pipe. The KURZ velocity probe was 
calibrated by the NASA Lewis Icing Research tunnel personnel prior to any engine inlet 
testing. 

The mass flows chosen for the Boeing 737-300 inlet tests were 23 lbm/s and 17.2 
lbm/s. These mass flows correspond to capture ratios of 1 and 0.75 respectively which are 
representative of flight conditions for this inlet. 

3.7 Droplet Impingement Data Acquisition 

To obtain the water droplet impingement data for a given model the following steps 
were executed during each impingement test. 


1. The pressure for the air and dyed water solution of the spray system were 
set to the required levels. These pressures were typically set to 106 psig for 
the dyed water at the supply tank (corresponding to 100 psig at the nozzle) 
and 65 or 80 psig for air. An air pressure of 65 psig produced a cloud with 
an MVD of 16 pim while the 80 psig air pressure produced a cloud with an 
MVD of 20 fim . 

2. The timer that controlled the duration of the spray was set to 2.6 seconds 
for spray clouds having an MVD of 20 pm and to 4.5 seconds for spray 
clouds with MVD of 16 pm . These time durations produced optimum dye 
impingement on the blotter strips for the concentration (0.0002 by weight) 
of the dyed water solution used in the 1989 tests. 

3. The blotter strips were attached at the required locations on the forward 
part of the test model. The blotter strips were typically 1.5 inches wide and 
of various shapes and lengths, depending on the model geometry and angle of 
attack. The strips were attached to the models with aluminum tape, which 
was pressed down smoothly to the model to minimize disturbances to flow 
in the locality of the blotter strips. The blotter paper used in the tests had 
a smooth and a rough side which yielded different optical properties. The 
smooth side was placed facing the flow. Each blotter strip was marked on 
its rough side to identify test run number, test model and strip location on 
the model. In addition, for the inlet tested the blotter strips were marked 
to indicate which part of the strip was attached to the outer surface and 
which part of the strip was attached to the inner surface of the inlet. Fig. 
3.7 shows typical of blotter strip installation on airfoils and inlet tested. 

4. The following test conditions were set prior to each impingement test: 

— spray system air and water pressures 

— spray duration 
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— angle of attack 

— temperature and humidity (approximately) 

— true air speed 

— mass flow (inlet only) was adjusted to the required level by 
opening or closing the regulating valve of the suction system. 

A precalibrated voltmeter was used to monitor inlet mass flow. 

5. The following test variables were recorded for each impingement test (see 
Appendix B): 

— dyed water concentration (constant for all tests) 

— test run and geometry tested 

— freestream static and total pressures (psia) 

— freestream average total temperature in °F measured at three 
different tunnel locations upstream of the spray bar system 

— true air speed 

— spray system water and air pressures (psig) 

— water pressure (psig) at the nozzles using electric transducers 

— spray time duration (seconds) 

— inlet mass flow in volts (converted to Ibm/s using mass flow 
calibration curve) 

6. The spray system was activated and a dye deposit was obtained on the blotter 
strips. The tunnel air speed was reduced to idle speed (about 10 mph) and 
the blotter strips were carefully removed from the model and placed in the 
control room and allowed to dry. The blotter strips were then placed in 
envelopes that were labeled for future reference. 

7. The model was wiped dean and new blotter strips were placed on the model 
for the next test. 

8. Steps 1 to 7 were repeated for all test runs performed. For each test condition 
five repeated runs were performed to obtain a statistical average. 

Test time between runs varied from 10 to 25 minutes depending on model and number 
of strips used. The tunnel was shut down completely between model changes. 
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FIGURE 3.1 

SCHEMATIC DIAGRAM OF SPRAY SYSTEM 



FIGURE 3.2 

NOZZLE ASSEMBLY DETAILS 
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FIGURE 3.3 

SPRAY SYSTEM INSTALLED IN IRT PLENUM CHAMBER 
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I -BEAM 



ALL DIMENSIONS IN iNCt-CS 


NOZZLE 
ASSEMBLY # 

NASA STANDARD 
NOZZLE # 

NOZZLE LOCATION 
COORDINATES IN INCHES 

Y 

z 

1 

40 

-17.5 

-0.5 

2 

86 

-29.25 

-29.25 

3 

27 

-5.0 

-31.25 

4 

109 

18.125 

-51.75 

5 

43 

37.25 

-29.375 

6 

50 

69.25 

-.25 

7 

48 

45.25 

-31.0 

8 

107 

19.75 

31.25 

9 

112 

-9.0 

36.0 

10 

111 

-30.5 

31.0 

11 

104 

23.25 

0.0 

12 

13 

-62.75 

-1.0 


FIGURE 3.4 

NOZZLE LOCATIONS WITH RESPECT TO IRT SPRAY BARS 


20 



(A) 1.25in x 1.25in BLOTTER SQUARES 



(B) 1.25in Wide x 3ft Long BLOTTER STRIPS 
FIGURE 3.5 

UNIFORMITY GRID WITH BLOTTER STRIPS ATTACHED 
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FIGURE 3.6 

REFERENCE COLLECTOR MECHANISM USED IN ALL 1989 TESTS 


(A) M5(l) 


0317 SWEPT AIRFOIL 



(B) 737 - 300 ENGINE INLET 
FIGURE 3.7 

BLOTTER STRIP LOCATIONS ON TEST MODELS 
(PAGE 1 OF 3) 
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(C) SWEPT WING TIP 
FIGURE 3.8 


BLOTTER STRIP LOCATIONS ON TEST MODELS 
(PAGE 2 OF 3) 
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(D) SMALL GLAZE ICE SHAPE 



(E) LARGE GLAZE ICE SHAPE 


FIGURE 3.8 

BLOTTER STRIP LOCATIONS ON TEST MODELS 
(PAGE 3 OF 3) 


25 


4.0 TEST MODELS AND CONDITIONS 


There were a total of six models tested. These included a Natural-Laminar-Flow airfoil 
NLF( 1) — 0414, a large and a small glaze simulated ice shapes, a swept medium speed 
M5(l) - 0317 airfoil, a swept wing tip with a NACA0012 airfoil section, and a 737 - 300 
engine inlet model. Descriptions of the test models and test conditions are given in this 
section and are summarized in Table 4.1. 

4.1 Two-Dimensional Flow Models 

4.1.1 N LF{ 1) — 0414 Airfoil Section 

This airfoil was designed in the early 1980’s and is representative of modern natural 
laminar flow airfoils. The NLF( 1) — 0414 was designed for 70% chord laminar flow on 
both surfaces, a lift coefficient of 0.4, a Reynolds number of 10 7 and a Mach number of 
0.4 (Ref. 13). 

The airfoil tested was made of mahogany with a fiberglass trailing edge, had a chord 
of 36 inches and spanned the height (6 ft) of the tunnel test section. Impingement data 
for this airfoil were obtained for two angles of attack 0 and 8 degrees and for two cloud 
MVD’s 16 and 20 microns. Two blotter strips A and B were attached to the airfoil during 
each test run. Strip A was located at midspan, 36 inches above the tunnel floor. Strip B 
wa s located 32 inches along the span from the tunnel floor. A summary of test conditions 
and number of repeats per condition is given in Table 4.1. The airfoil geometry definition 
is given in Figs. 4.1a and 4.1b. Figure 4.2a shows the installation of this airfoil in the IRT 
test section. 

4.1.2 Simulated Ice Shapes 

The large and small glaze simulated ice shapes tested are representative of the types 
often encountered on aircraft surfaces. Each ice shape was 18 inches in span and was 
mounted on a two inch in diameter and 6 ft long cylinder as shown in Fig. 4.2b. The only 
test parameter varied for the impingement tests of the two ice shapes was cloud MVD. 
Test conditions for these geometries are presented in Table 4.1. The coordinates of the ice 
shapes are given in Figs. 4.1c and 4. Id. 

4.2 Three-Dimensional Flow Models 

Models with three-dimensional geometry and/or three-dimensional flow fields tested 
in 1989 included an infinite swept wing, a low aspect ratio swept wing and an engine inlet. 
All geometries tested are representative of modern aircraft surfaces. 
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4 . 2.1 MS( 1) —0317 Swept Airfoil 

The MS(1) — 0317 wing section is representative of medium speed airfoils. It was 
designed in the mid 1970s for general aviation aircraft (Ref. 14). The airfoil was designed 
for a lift coefficient of 0.3 a Reynolds number of 1.4xl0 7 and a Mach number of 0.68. 

The model tested was 6 ft span 30 degree swept wing made of mahogany with a 36 
inch chord in the streamwise direction. The coordinates of the AfS(l) — 0317 section are 
given in Fig. 4.3a. Model installation in the tunnel test section is shown in Fig. 4.4a. 

Impingement data for this geometry were obtained for two angles of attack, 0 and 8 
degrees, and two MVD’s 16 and 20 (im . Blotter strips were attached to the airfoil surface 
at locations A and B, 36 and 32 inches above the tunnel floor respectively. Each test 
condition was repeated five times. 

A 36 inch chord unswept MS(l) - 0317 airfoil with 6 ft span was tested in 1985. 
The swept model used in the 1989 tests was chosen to investigate how sweep angle affects 
impingement characteristics. 

4 . 2.2 Swept Wing Tip 

This was a wooden finite wing model with 30° leading edge sweep. The model had 
a NACA 0012 wing section normal to the leading edge with a chord length of 15 inches. 
The wing span was 15 inches and the aspect ratio of the wing was 1.4. This model was 
selected to study how droplet impingement is affected by the wing tip vortex. 

Impingement tests for this model were conducted at two angles of attack, 0 and 8 
degrees, and for two MVD sizes, 16 and 20 /xm . Blotter strips were attached to this 
model at three locations. Location A was at the tip of the wing, location B was 42 inches 
above the tunnel floor and location C was 32 inches above the tunnel floor. All strips were 
attached to the leading edge of the wing and were shaped to follow the wing surface along 
the streamwise direction. Blotter strips B and C were V-shaped to fit the wing leading 
edge. Installation of the wing tip model with blotter strips attached in the IRT test section 
is shown in Fig. 4.4b. 

4 . 2.3 Boeing 737-300 Engine Inlet Model 

This 0.2547 scale inlet model (Boeing Model 5718-1-1) was made out of solid aluminum 
and was provided by the Boeing Company. Details of this geometry are available on 
magnetic tape from NASA Lewis Research Center. A brief geometry definition for this 
inlet is given in Fig. 4.3. Inlet areas of interest are: fan face area 149.77 in 2 , throat area 
143.01 in 2 and highlight area 185.58 in 2 . 

The 737-300 inlet was tested at inlet mass flows of 17 and 23 Ibm/sec (corresponding 
to capture area ratios of 0.75 and 1 respectively), angles of attack of 0 and 15 degrees, 
and MVD’s of 16 and 20 /xm , for a total of four test conditions (see Table 4.1). Each 
test condition was repeated several times. Blotter strips were attached to the inlet at 
intervals of 45 degrees giving a total of eight strips. The maximum angle of attack was 
limited to 15 degrees to keep the inlet model within the uniform cloud region. The inlet 
was mounted in the tunnel section so that the lower engine lip was parallel to the tunnel 
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left wall looking downstream. Variations in angle of attack were obtained by rotating the 
inlet about an axis perpendicular to the tunnel turn table. The inlet installation with 
mass flow connection and inlet support is shown in Fig. 4.4b. Blotter strips attached to 
the inlet surface are also shown in this figure. 

Impingement data for the 737-300 inlet were also obtained in 1985. The inlet was 
retested in 1989 in order to resolve discrepancies, between the 1985 experimental and 
analysis data, observed at a = 15° near the inlet upper lip. These discrepancies were 
attributed to cloud non-uniformity. 

4.3 Reference Collector Mechanism 

The reference collector was tested at two tunnel locations corresponding to the location 
of the engine inlet at 0 and 15 degrees angle of attack. Two droplet MVD’s, 16 and 20 fim 
were used in the collector impingement tests. Each collector test condition was repeated 
at least five times. Nine blotter strips were attached to the collector as shown in Fig. 4.4d. 
The long center blade of the collector was used to provide LWC information for all wings 
and ice shapes. This blade was used in all collector tests with the collector at the a = 0° 
location. The eight shorter blades were used to determine LWC for the engine inlet. 
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TABLE 4.1 

TEST CONDITION MATRIX 


MODEL DESCRIPTION 

TAS 

mph 

a 

(deg) 

MVD 

pm 

Mass Flow 
lbm/s 

Blotter 

Strips Per Run 

Test Runs 
(See App. B) 

Simulated Iced Shapes 

165 

- 

20 

— 

2 

116- 126 

NLF(1)-0414F Airfoil, 36 inch Chord 

165 

0 , 8 

16 , 20 

- 

2 

215 - 236 

30° Swept MS(1)-0317 Airfoil, 36 inch 
Chord 

165 

0 , 8 

16 , 20 

- 

2 

148 - 167 

30° Swept NACA 0012 Wing Tip, 15 
inch Chord 

165 

0 , 8 

16 , 20 

- 

3 

127 - 147 

Boeing 737-300 Inlet Model 

171 

0 , 15 

16 

17 

8 

168 - 173, 179- 180 

Boeing 737-300 Inlet Model 

173 

0 

20 

23 

8 

174 - 176, 177 - 178 

Boeing 737-300 Inlet Model 

171 

15 

20 

23 

8 

174 - 176, 177 - 178 

Reference Collector Mechanism 

165 

0 , 15 

16 , 20 

- 

8 - 9 

95 - 115 


x/c 

Y/C 


X/C 

Y/C 


x/c 

Y/C 

. 000000 

. 000000 

.245187 

.074349 

.735392 

.047492 

. 000085 

. 001585 

.261054 

.075830 

.750058 

. 042542 

. 000299 

. 003274 

.277233 

.077161 

.764925 

. 037208 

. 001231 

. 007144 


.293699 

.078380 


.779951 

. 031694 

. 002695 

. 010618 

.310424 

.079454 

.795034 

. 026178 

. 004989 

. 014163 

.327391 

. 080369 


.810124 

. 020750 

. 008005 

.017552 

.344571 

.081151 

.825179 

.015483 

.011774 

. 020769 

.361925 

.081781 

.840076 

. 010464 

. 016268 

. 023816 

. 379421 

.082240 

.854693 

.005783 

. 021468 

. 026795 

.397052 

. 082536 

.868960 

. 001467 

. 027356 

.029735 

.414812 

.082677 

.882768 

-.002475 

. 033891 

. 032633 

.432667 

.082633 

.896006 

-.006044 

. 041042 

. 035480 

.450558 

.082429 

.908644 

-.009267 

. 048811 

. 038317 

.468450 

.082047 


.920659 

-.012161 

.057201 

. 041092 

.486327 

.081507 

.931980 

-.014739 

. 066189 

. 043825 

. 504159 

.080794 


.942511 

-.017008 

L 07 57 67 

. 046482 

. 521931 

.079893 

.952200 

-.018994 

. 085915 

. 049070 

.539641 

.078779 

.961042 

-.020722 

.096610 

.051588 

.557254 

.077489 

.969034 

-.022206 

. 107826 

.054033 

. 574742 

.075988 


.976155 

-.023456 

. 119545 

. 056398 


.592064 

.074285 


.982370 

-.024492 

. 131756 

.058692 

.609177 

.072377 


.987660 

-.025333 

. 144443 

.060917 

.626040 

.070245 


.992021 

-.026006 

. 157592 

. 063092 

.642629 

.067900 


.995456 

-.026519 

. 171193 

.065206 

.658928 

. 065348 


.997952 

-.026872 

. 185212 

.067240 

.674926 

.062510 


.999480 

-.027067 

. 199628 

.069172 

.690586 

.059376 


1.000000 

-.027122 

.214447 

.071009 

.705860 

.055889 



.229647 

.072735 

.720751 

.055194 




(A) NLF{1) - 0414 AIRFOIL (UPPER SURFACE) 
FIGURE 4.1 

DEFINITION OF TWO-DIMENSIONAL TEST GEOMETRIES 

(PAGE 1 OF 4) 
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x/c 

Y/C 

.000000 

.000000 

. 000085 

-.001535 

.000164 

-.002120 

. 000740 

-.004536 

.002095 

-.006984 

. 004175 

-.009008 

.007129 

-.010993 

. 010874 

-.012933 

.015540 

-.014882 

.021096 

-.016854 

.027380 

-.018787 

. 034569 

-.020742 

. 042393 

-.022654 

.050985 

-.024572 

.060274 

-.026487 

.070243 

-.028383 

.080881 

-.030259 

.092159 

-.032116 

. 104058 

-.033945 

. 116557 

-.035741 

. 129635 

-.037497 

. 143277 

-.039212 

. 157457 

-.040888 

. 172148 

-.042421 

. 187328 

-.044107 

.202969 

-.045646 

.219043 

-.047125 


x/c 

Y/C 

. 726433 

-.046065 

.743743 

-.042296 

.761642 

-.038850 

.779550 

-.035991 

.797188 

-.033529 

.814513 

-.031444 

.831368 

-.029735 

.847719 

-.028310 

. 863493 

-.027230 

. 878523 

-.026450 | 

. 892802 

-.025925 

.906336 

-.025641 | 

.919043 

-.025539 | 

.930841 

-.025569 | 

.941715 

-.025689 

.951668 

-.025861 

.960696 

-.026061 

.968804 

-.026275 

.975996 

-.026483 | 

.982266 

-.026675 

.987613 

-.026858 

.992033 

-.027036 

. 995503 

-.027211 

.997994 

-.027367 

.999497 

-.027475 

1.000000 

-.027514 




x/c 

Y/C 

.235525 

-.048542 

.252387 

-.049901 

.269586 

-.051189 

.287087 

-.052411 

.304866 

-.053561 

. 322901 

-.054635 

.341156 

-.055635 

.359611 

-.056539 

. 378260 

-.057344 

.397074 

-.058052 

.416017 

-.058658 

.435049 

-.059142 

.454127 

-.059517 

.473222 

-.059785 

.492319 

-.059950 

.511402 

-.060012 

.530430 

-.059979 

.549361 

-.059792 

.568160 

-.059456 

.586782 

-.058982 

. 605204 

-.058340 

. 623397 

-.057533 

. 641303 

-.056524 

. 658920 

-.055246 

. 676262 

-.053698 

.693229 

-.051845 

.709795 

-.049388 


(B) NLF{ 1) -0414 AIRFOIL (LOWER SURFACE) 


FIGURE 4.1 

DEFINITION OF TWO-DIMENSIONAL TEST GEOMETRIES 

(PAGE 2 OF 4) 
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COORDINATES FOR 2 MINUTE GLAZE ICE ACCRETION 
(NOTE : DATA FROM NASA - LEWIS) 


C = .0508 m = 2.0" 





X /c 





V 

1.0000000 

0.9980649 

0.9922500 

0.9825492 

0.9689901 

0.9516240 

0.9310589 

0.9071791 

0.8801693 

0.8502499 

0.8188995 

0.7850649 

0 .7489684 

0.7108700 

0.6710334 

0.6297538 

0.5873248 

0.5440689 

0.^999999 

0.4759901 

0.4524940 

0.4290000 

0.4054999 

0 . 3820097 

0 . 3584803 

0.3355491 

0.3122243 

0.2890708 

0.2673995 

0.2446043 

0.2220708 

0.1991948 

0.1762047 

0.1531397 

0.1303445 

0 .1108405 

0 . 0 909448 

0.0724606 

0 .0558996 

0.0394547 

0.0247559 

0 .0099842 

- 0.0042106 

- 0.0170197 

- 0.0281004 

- 0.0373091 

-0 . 0450689 

- 0 . 0515197 

- 0 .05801 97 

- 0 . 0644389 

- 0.0698149 

- 0.0747105 

- 0.0804842 

-0 . 0849507 

-0 . 0852500 

- 0 . 0850000 

- 0 . 0850000 

- 0.0854055 

- 0.0842637 

- 0 . 0799999 

- 0.0839192 

- 0 . 0852854 

- 0 . 0850196 

- 0.0853640 

- 0.0842637 

-0 . 081 3838 

- 0.0760491 

- d .0694646 

-0 . 0637953 

- 0.0585098 

- 0 .0538406 

-0 . 0481850 

-0 . 0410039 

- 0.0333346 

- 0.0256988 

- 0.0163307 

- 0.0048189 

0 .0073898 

0.0204154 

0.0348799 

0 .0477657 

0.0610098 

0.0794842 

0.0971141 

0.1165746 

0 .1383858 

0.1608444 

0.1835490 

0.2056909 

0.2272854 

0.2490747 

0.2716141 

0.2944901 

0.3174704 

0.3403739 

0.3624055 

0.3854153 

0.4083858 

0.4313503 

0.4542657 

0.4771495 

0.4999999 

0.5395039 

0.5788838 

0.6178936 

0.6562834 

0.6938090 

0.7302244 

0.7645295 

0.7971752 

0.8279586 

0.8566948 

0.8823385 

0.9058799 

0 . 9271240 

0.9459153 

0.9620984 

0.9755492 

0.9861495 

0.9938051 

0 .9984390 

1.0000000 

0 .00 000 00 

0.0000000 

0.0000000 

0 .0 000 000 


Y /c 






0 . 0000000 

- 0.0439291 

- 0.0876947 

- 0.1309448 

-0 . 1733444 

- 0.2145551 

- 0.2533504 

- 0.2901789 

- 0.3247655 

- 0.3568247 

- 0.3851002 

- 0.4107795 

- 0.4336042 

- 0.4533602 

- 0.4698345 

- 0.4528699 

- 0.4923149 

- 0.4980550 

- 0.4999999 

- 0.4999999 

- 0.4999999 

- 0.4999999 

- 0.4999999 

- 0.4999999 

- 0.5004880 

- 0.4961358 

- 0.4949743 

- 0.4965157 

- 0.5057697 

- 0.5114841 

- 0.5181849 

- 0.5234605 

- 0.5227931 

- 0.5178602 

- 0.5121495 

- 0.4995293 

- 0 .4870098 

- 0.4728503 

- 0.4561396 

- 0.4393857 

- 0.4210649 

- 0.4027755 

- 0 . 3840905 

- 0.3643640 

- 0.3437046 
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(C) SMALL GLAZE ICE SHAPE 
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COORDINATES FOR 15 MINUTE GLAZE ICE ACCRETION 
(NOTE DATA FROM NASA - LEWIS) 
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2 0 ” 


X /r 
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0.99558*5 

0 . 9822*99 

0.9605*52 

0.9307696 

0.893*802 

0.8*93739 

0.8019350 

0.7*90787 

0.69157*8 

0.63032*7 

0.5662991 

0 .*999999 

0. *750000 

0. **99999 

0 .*250000 

0.3999999 

0.37*9999 

0.3515000 

0.3279999 

0 . 30**999 

0.2810000 

0.257*999 

0.2339998 

0.210*999 

0.1869999 

0.1635000 

0.1399999 

0.116*998 

0.0929999 

0 . 0695000 
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-0 .1890000 
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-0 .*00*999 

-0.4239998 
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-0.4599999 

-0. *800000 

rO . *899999 

-0.4899999 

-0.4899999 

-0.4750000 

-0.4899999 

-0 .*899999 

-0 .'*899999 

-0.4800000 
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-0.4800000 
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-0.49**999 

-0.4710000 
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-0.4335*51 
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-0.555*978 
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(D) LARGE GLAZE ICE SHAPE 
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(A) NLF{ 1) - 0414 AIRFOIL 
FIGURE 4.2 


INSTALLATION OF TWO-DIMENSIONAL MODELS IN IRT TEST SECTION 

(PAGE 1 OF 2) 


34 



(B) SMALL GLAZE ICE SHAPE 



(C) LARGE GLAZE ICE SHAPE 


FIGURE 4.2 

INSTALLATION OF TWO-DIMENSIONAL MODELS IN IRT TEST SECTION 

(PAGE 2 OF 2) 
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x/c 

Y/C 

UPPER 

Y/C 

LOWER 

.50000 

. 09535 

-.06591 

.52500 

. 09323 

-.06389 

. 55000 

. 09073 

-.06138 

.57500 

. 08777 

-.05845 

. 60000 

. 08448 

-.05501 

. 62500 

. 08079 

-.05106 

. 65000 

. 07672 

-.04674 

.67500 

. 07232 

-.04214 

.70000 

. 06763 

-.03735 

. 72500 

. 06269 

-.03255 

.75000 

. 05755 

-.02780 

.77500 

. 05225 

-. 02309 

.80000 

. 04687 

- . 01857 

. 82500 

. 04132 

-.01433 

. 85000 

. 03576 

-.01049 

. 87500 

. 03013 

-.00719 

.90000 

. 02444 

-.00460 

.92500 

.01873 

-.00289 

.95000 

. 01302 

-.00232 

.97500 

. 00720 

-.00324 

1.00000 

. 00125 

-.00597 







1 : 




x/c 

Y/C 

UPPER 

Y/C 

LOWER 

. 00000 

. 00099 

.00099 

. 00200 

.01248 

-.00857 

. 00500 

.01950 

-.01366 

. 01250 

. 03099 

-.02105 

. 02500 

. 04322 

-.02866 

. 03750 

. 05210 

-.03423 

. 05000 

. 05893 

-.03865 

. 07500 

.06840 

-.04541 

. 10000 

. 07511 

-.05058 

. 12500 

. 08033 

-.05477 

. 15000 

. 08454 

-.05817 

. 17500 

. 08805 

-.06099 

.20000 

. 09096 

-.06330 

. 22500 

.09339 

-.06527 

. 25000 

.09536 

-.06685 

. 27500 

. 09694 

-.06812 

• 30000 

. 09815 

-.06909 

.32500 

.09901 

-.06978 

.35000 

. 09952 

-.07021 

.37500 

.09972 

-.07036 

.40000 

.09956 

-.07019 

.42500 

.09909 

-.06967 

.45000 

.09826 

-.06880 

| . 4 75C0 

.09700 

- 06755 


(A) MS{ 1) -0317 30° SWEPT AIRFOIL 
FIGURE 4.3 

DEFINITION OF THREE-DIMENSIONAL TEST GEOMETRIES 

(PAGE 1 OF 5) 
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X 

(per cent c) 

V 

(per cent c) 

0 

0 

0.5 


1.25 

1.894 

2.5 

2.615 

5.0 

3.555 

7.5 

4.200 

10 

4.683 

15 

5.345 

20 

5.737 

25 

5.941 

30 

6.002 

40 

5.803 

50 

5.294 

60 

4.563 

70 

3.664 

80 

2.623 

90 

1.448 

95 

0.807 

100 

0.126 




(B) NACA 0012 30° SWEPT WING TIP 
FIGURE 4.3 

DEFINITION OF THREE-DIMENSIONAL TEST GEOMETRIES 

(PAGE 2 OF 5) 
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(C) 737-300 BOEING ENGINE INLET - FRONT VIEW 

FIGURE 4.3 

DEFINITION OF THREE-DIMENSIONAL TEST GEOMETRIES 

(PAGE 3 OF 5) 
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All dimensions in inches 

(D) 737-300 BOEING ENGINE INLET - CROSS SECTION 

FIGURE 4.3 

DEFINITION OF THREE-DIMENSIONAL TEST GEOMETRIES 

(PAGE 4 OF 5) 
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S ec-A A 

(E) 737-300 BOEING ENGINE INLET - CENTER BODY 

FIGURE 4.3 

DEFINITION OF THREE-DIMENSIONAL TEST GEOMETRIES 

(PAGE 5 OF 5) 


(A) MS(1) - 0317 30° SWEPT AIRFOIL 



(B) NACA 0012 30° SWEPT WING TIP 


FIGURE 4.4 


INSTALLATION OF THREE-DIMENSIONAL TEST MODELS 
IN IRT TEST SECTION 
(PAGE 1 OF 2) 
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(C) 737-300 BOEING ENGINE INLET 
FIGURE 4.4 


INSTALLATION OF THREE-DIMENSIONAL TEST MODELS 
IN IRT TEST SECTION 
(PAGE 2 OF 2) 
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FIGURE 4.5 

AXIAL LOCATIONS OF TEST MODELS IN IRT TEST SECTION 






5.0 DATA REDUCTION METHOD 


The laser reflectance method used to reduce the 1989 impingement data is detailed 
in Ref. 10. In summary, the reflectance method is based on the assumption that the 
reflectance of the dyed paper is a measure of the dye mass per unit area of the paper. 
Dyed regions on the blotter strip corresponding to high impingement rates are darker in 
color and reflect less light than those corresponding to low impingement rates. Regions 
with no dye accumulation are white and reflect the maximum amount of laser light. The 
relation between dye concentration and laser reflectance is not linear and is determined 
from calibration tests. To maximize the sensitivity of the reflectance method, the dye used 
in the tests must have a strong absorption at the wavelength of the laser radiation. This is 
accomplished by using a blue dye whose maximum absorption occurs at 629.5 nm, which 
is very close to the wavelength (632.8 nm) of the red He-Ne laser employed in the data 
reduction system. 

5.1 Reflectance Calibration Curve 

The laser reflectance calibration curve relates laser reflectance to dye mass on the 
blotter strip. The curve is a standard against which the reflectance of each blotter strip is 
compared during the data reduction process. If the type of dye, the blotter paper and the 
laser instrumentation remain constant the curve does not change. 

Due to minor changes made to the laser reflectometer system a new calibration curve 
was produced for the analysis of all 1989 test data. To produce the new laser reflectance 
calibration curve a series of tests were performed in the IRT during the 1989 test period. 
In each of these tests a large square of blotter paper was placed on the uniformity grid 
and was exposed to the spray cloud for a short period of time. By varying the time 
exposure to the spray cloud blotter squares with various color intensities were obtained. 
Blotter samples were cut from each large blotter square. The reflectance and dye mass for 
each of the samples was measured using the laser and absorption spectroscopy methods 
respectively. The dye mass from each blotter sample was then divided by the area of the 
sample and the reflectance of each sample was normalized by the reflectance of the white 
blotter paper. A total of 19 dyed blotter samples of various areas were used to define the 
1989 laser calibration curve. Raw and normalized calibration data are given in table 5.1. 
The 1989 and 1985 laser calibration curves are shown in Fig. 5.1. Both calibration curves 
show similar trends. 

5.2 Automated Data Reduction System 

The laser reflectometer and data acquisition system used for the analysis of the 1985 
impingement data were also utilized for the current tests. A block diagram of the complete 
data reduction system is presented in Fig. 5.2. Details of the design aspects of the system 
are given in Ref. 10. A brief description of the data reduction system is provided below. 

The main components of the reflectometer are (a) a red He-Ne laser used to scan the 
blotter strips, (b) a rotating drum for mounting the blotter strips, (c) a convergent lens 
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for collecting the reflected light, and (d) a photodetector which converts the reflected light 
collected by the lens into a voltage (Vj) which is stored for further analysis. A splitter 
glass plate and another photodetector are also shown and are used to monitor fluctuations 
in laser output. The voltage (Vj) from the second photodetector is also stored and is used 
in the data analysis. 

The digital data acquisition system controls the operation of the reflectometer and is 
used to analyze and plot the data. The main components of this system are a desk top 
computer with two disk drives, a plotter and a printer, and two voltmeters for monitoring 
the voltages from the photodetectors. 

5.3 Data Reduction Process 

The reduction of the impingement data was performed in two main steps. During the 
first step reflectivity values were extracted from the dyed blotter strip and were stored in the 
form of voltages on disk. During the second step the reflectivity values were transformed 
into plots of impingement efficiency versus surface distance measured with respect to some 
reference point on the model. 

To extract the impingement data from a single blotter strip, the strip was mounted on 
the drum of the laser reflectometer and one to three horizontal scans were obtained depend- 
ing on blotter strip geometry. Straight strips were scanned at three different horizontal 
locations. Curved strips were scanned at one horizontal location. Curved strips were re- 
quired for the swept MS-0317 airfoil and the swept wing tip in order to fit these geometries 
along constant streamwise sections. Because the whole strip could not be mounted on the 
drum these strips had to be scanned in small sections and hence were much more difficult 
to analyze. The data from each section of each curved strip were stored in different files 
which had to be combined into a single file for final analysis (see Appendix A). 

For each scan the voltages V\ and Vj* from the two photodetectors were recorded and 
stored on disk. The reflectometer takes 47 reflectivity readings per inch of each horizontal 
scan at equal space intervals. 

To transform reflectivity into a plot of impingement efficiency versus surface distance, 
the following steps were used: 


1. The raw reflectivity values were divided by the reflectivity of the bare blotter 
paper to obtain normalized reflectivity data. The following equation was used 
for each data point of the dyed blotter paper 


Jotter Paper 
( Vl /F 2 ) BareBlotterPaper 


Note that {VifV 2 ) BareBlotterPaper is the average of 1000 data points mea- 
sured on a reference blotter strip. 

2. The normalized reflectivity data were converted into dye mass per unit area 
using the laser calibration curve. 
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3. The impingement efficiency value for each data point recorded was obtained 
from the following equation. 


0 


Dye mass per unit area corresponding to a given data point 


X 0 Collector 


Average mass per unit area of reference collector 

(5 - 2) 

This equation is valid if the collector and the model have been exposed to 
sprays of identical duration, as is the case here. 


All collector strips were analyzed prior to the analysis of the test model strips. Each 
collector strip was scanned and an average mass per unit area was calculated for the strip. 
This value was used as the denominator in eq. (5-2). The value of 0coiiector depends on 
the MVD of the spray cloud and on cloud droplet distribution. The MVD’s of the spray 
clouds used in the current tests were 16 pm and 20pm. The corresponding theoretical 
collector efficiency values were 0.86 and 0.89, 


5.4 Laser Reflectance Method Error Analysis 

The repeatability and accuracy of the laser reflectance method depends on the accu- 
racy of the laser reflectance calibration curve, the stability of the laser output power, the 
uniformity of blotter paper composition, the degree of dye penetration into the blotter 
paper, and the dye concentration range over which data are collected. These factors are 
discussed below. 

The dye concentration values for the data points used in defining the laser reflectance 
calibration curve were obtained by absorption spectroscopy. The accuracy of this method 
was verified by analyzing several samples of dye solution of known concentrations. These 
concentrations were in the range of 200 to 500 milligrams of dye per liter of water. Ab- 
sorption spectroscopy reproduced the concentrations of the sample solutions accurate to 
within .6 mllligram/liter (i.e., ±0.3%). 

To further validate the calibration method, groups of dyed blotter samples with nearly 
identical reflectance values were produced and analyzed by absorption spectroscopy. Re- 
sults from this test are shown in Table 5.1. 

The relation between normalized reflectance and dye mass is monotonic; dye mass per 
unit area increases as reflectance is decreased. The results of Table 5.1 can be used to 
determine an estimate for the uncertainty in the calibration curve. Upon close examina- 
tion of the data it is seen that in some cases blotter samples with identical reflectivity 
values differ in dye mass by less than 2% from the average mass. Also, samples for which 
absorption spectroscopy produced identical dye concentrations differ in reflectivity by less 
than 2% from the average reflectivity. It is reasonable to conclude that the uncertainty in 
the calibration curve is of the order of 2%. 

During the early stages of the development of the laser reduction method, it was found 
that the output light intensity of the laser beam varied significantly over long periods of 
time. Since the light reflected from the blotter paper is a direct function of the incident 
light, the laser fluctuations must be monitored and taken into account during the data 
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reduction. First, the He-Ne laser used in the preliminary studies was replaced with a 
polarized He-Ne laser. Experience indicates that polarized lasers give a more stable out- 
put. Second, optical feedback was incorporated to monitor any remaining fluctuations in 
laser power. During the data reduction, the effect of laser power fluctuations is essentially 
eliminated by using the ratio of the collected signal to the feedback signal as explained 
in section III. The optical feedback design was tested by scanning a reference bare blot- 
ter strip a number of times over a period of several months. Each time 400 reflectance 
measurements were obtained and the average of the 400 data points was evaluated. The 
standard deviation of these averaged reflectance values was found to be ±0.3%. 

The uniformity of the reflectivity of an unsprayed (bare) piece of blotter paper was 
measured by scanning 200 points at five horizontal locations. The variation in reflectance 
was ±1.3% from the average value based on 1000 (5 x 200) reflectance measurements. This 
variation in reflectance is due to the roughness of the paper itself and due to variation in 
the composition of the blotter paper. 

Penetration of dye into the blotter paper should be avoided because the data reduction 
is based on surface reflectance measurements. This was accomplished by keeping the spray 
time short. A series of tests were performed to establish the optimum spray time for 
each droplet size used in the actual impingement experiments. During these tests blotter 
strips were exposed to sprays of various time durations. After each test the penetration of 
dye into the paper was examined using a microscope. Spray times of the order of 3 to 5 
seconds produced penetrations of the order of 10% of the blotter thickness. This level of 
penetration was acceptable and did not affect the quality of the data. 

The discussion above indicates that the error in measuring reflectance is of the order 
of 2 %. The corresponding error in extracting the dye mass per unit area of blotter paper 
needs to be determined. Let the absolute error in reflectance be dR . The absolute error in 
dye mass per unit area dC is a function of R and dR where R is the normalized reflectance. 
For R values close to 1 a given dR has a small dC associated with it. Near the low end 
of the reflectivity axis (e.g. for R values near 0.3) the same value of dR produces a large 
dC. In both cases the percent error in dye mass per unit area {dCjC) is large. In the first 
case dC is small, but C is also very small, so the ratio dC/C is large. In the second case 
dC is large, while C has some finite size, so the ratio dC fC remains large. Between these 
extremes there is a region where dCjC is reduced to a minimum. This region of low values 
in dCjC corresponds to a certain range of absolute normalized reflectivity. This range 
of reflectivity is determined by differentiating R with respect to C and plotting dC/C as 
a function of R. The calibration curve is first approximated with cubic splines and the 
differentiation is performed numerically. Plots of dC(C versus R for dR = 0.01, 0.02 and 
0.03 are given in Fig. 5.3. To maintain, for example, the percent error in the measurement 
of dye mass per unit area to less than 10%, assuming 1% error in normalized reflectance 
(i.e., dR = 0.01), the normalized reflectance of the dyed blotter paper must be in the range 
of 0.57 to 0.90. We also note that to obtain reflectance values within the required range the 
concentration of the dye solution and the duration of the spray must be chosen carefully. 


48 



5th Degree Least Squares Polynomial □ Normalized Raw Data 
R = - 071931m 5 +0,5398m 4 -L5759m 3 +2.2269m 2 -i.5934m+1.0i09 
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STANDARD LASER CALIBRATION CURVE 
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FIGURE 5.2 

BLOCK DIAGRAM OF DATA REDUCTION SYSTEM 



FIGURE 5.3 

PLOT OF dCjC AS A FUNCTION OF NORMALIZED REFLECTIVITY 
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6.0 ANALYSIS METHODS 

Computer codes to predict the local impingement efficiency of two-dimensional geome- 
tries have been available since the 1940s. More recently, three-dimensional codes have been 
developed to predict the local impingement efficiency distribution of complex geometries. 
These codes are currently used in the design of de-icing and anti-icing systems of modern 
aircraft. 

Analysis data presented in this report have been obtained by two- and three- 
dimensional trajectory computer codes. In these codes the calculation of water impinge- 
ment involves the following main steps: 


• Body surface definition 

• Calculation of the flow field 

• Droplet trajectory calculation 

• Calculation of local impingement efficiency as a function of surface distance. 

Brief descriptions of the computer codes used to produce the analysis impingement 
data presented in Section 7 are given below. 

6.1 Two-Dimensional Analysis Code 

A trajectory computer code from Boeing Military Airplanes (Ref. 15) were used to 
produce analysis data for the NLF-0414 airfoil. In this code the surface geometry definition 
is a cubic spline fit between user specified points. The flow field is generated using a finite 
difference formulation of the compressible streamfunction equation. Droplet trajectories 
are calculated by solving equation (2 - 1) by a modified Euler scheme, for each droplet size 
of the discrete distribution given in Fig. 2.1. Local impingement efficiency is calculated 
using the a two-dimensional version of the method described in Section 6.2 below. 

6.2 Three-Dimensional Analysis Code 

Analysis impingement data for the Boeing 737-300 inlet were obtained using a Boeing 
trajectory computer code. The three-dimensional analysis follows the same basic steps 
as the two-dimensional analysis code described in Section 6.1.1. The aircraft surface is 
defined by bi-cubic splines that provide a smooth and well defined surface (Ref. 16). The 
flow field is obtained by a full potential code (Ref. 17). The potential flow velocity is 
calculated at both mesh intersections and mesh-surface intersections. 

Individual water droplets, whose size is determined from the distributions measured 
experimentally (see Fig. 2.1), are traced to their point of impingement on the body surface. 
Droplet trajectories are calculated by solving equation (2 - 1) by a variable step 4th order 
Adams predictor corrector method (Refs. 1,18). At each computed trajectory position, 
the potential flow computational mesh cell is located and the flow velocity at that position 
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is linearly interpolated from the cell corner velocities. This is true for all flow field cells. 
For cells intersecting the body surface a least squares technique is used to determine local 
velocities. 

The impingement efficiency j3 is evaluated numerically as the ratio of impinging surface 
particle flux to the freestream flux. The local value of f3 for each point on the body surface is 
simply the ratio of that infinitesimal area at freestream to the corresponding impingement 
area at the surface as shown in Fig. 2.2. 
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7.0 RESULTS AND DISCUSSION 


Experimental water droplet impingement data for all geometries tested are presented in 
this section. Analysis impingement data are also presented for the NLF(1)-0414 airfoil and 
the Boeing 737-300 engine inlet. All data are in the form of local impingement efficiency 
plotted versus surface distance 0 - curve). Surface distance is measured with respect to a 
reference point termed the highlight. For all test models except ice shapes, the highlight 
point is the most forward (upstream) point on the geometry at an angle of attack of zero 
degrees. This point is also used as the reference point for all other angles of attack. For 
airfoil and wing geometries the surface distance is negative along the upper surface and 
positive along the lower surface. For the engine inlet, positive surface distance is along the 
inlet inner surface while negative surface distance corresponds to the inlet outer surface. 
The highlight locations and surface distance conventions for airfoils, ice shapes and inlet 
are summarized in Figs. 7.1 and 7.2. The main features of each - curve are: 

1. Curve overall shape 

2. Magnitude and location of peak efficiency 

3. Extent of impingement limits 

AH experimental curves presented were obtained by averaging data from several blotter 
strips. In some cases experimental data from as many as ten blotter strips were used to 
produce a single averaged impingement curve. Consider for example the large glaze ice 
shape. Five test runs were performed for this geometry at identical test conditions. During 
each run two blotter strips were attached to the surface of the ice shape giving a total 
of ten blotter strips per test condition. For most geometries five test runs at identical 
test conditions were performed (see Table 4.1). By repeating each test several times and 
averaging the results, small fluctuations in tunnel conditions and spray system performance 
are averaged out. 

The blotter strips used to produce each averaged impingement curve are given at the 
top of each curve. The notation used to identify each strip is as follows: 

RUN G E AA BS 


where, 

1. RUN indicates the run number 

2. G indicates the type of geometry. Possible values for G are: N for NLF{ 1) — 
0414, 1 for Ice shapes, M for MS(1)-0317, W for Wing Tip and 3 for 737-300 
inlet. 

3. E indicates the collector efficiency. A value of 6 corresponds to a collector 
efficiency of 0.86 and an MVD of 16 fim . A value of 9 corresponds to a 
collector efficiency of 0.89 and an MVD of 20 fim , 
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4. AA is the Angle of Attack with possible values of 0, 8 and 15 degrees. 

5. BS indicates the blotter strips used in the averaging process. For each geom- 
etry tested blotter strips were attached to the geometry at several locations 
identified by capital letters. Blotter strips used and corresponding locations 
for all geometries tested are given below. 

a . NLF( 1) — 0414 - Blotter strips A and B corresponding to lo- 
cations along the airfoil span at 36 and 32 inches above the 
tunnel floor respectively. 

b . ICE SHAPES - Blotter strips A and B corresponding to lo- 
cations along the span of each ice shape at 36 and 32 inches 
above the tunnel floor respectively. 

c . M5(l) — 0317 - Blotter strips A and B corresponding to loca- 
tions along the wing span at 36 and 32 inches above the tunnel 
floor respectively. 

d . WING TIP - Blotter strip A at the wing tip and blotter strips 
B and C corresponding to locations along the wing span at 42 
and 32 inches above the tunnel floor respectively. 

e . 737-300 engine inlet - Blotter strips A, B, C, D, E, F, G and 
H at 45 degrees intervals from the inlet upper lip. Strip A 
corresponds to 6 = 0° and strip H corresponds to 6 ~ 315° 

(see Fig. 3.7B). 


Consider for example Figure 7.29A2. The strip identification numbers shown at the 
top of this figure are 1773900BH and 1783900BH. These IDs indicate that four strips were 
used to produce the average impingement curve, the B and H strips of runs 177 and 178 
of the 737-300 inlet tested at an angle of attack of 0 degrees and an MVD of 20 /im . 

All experimental data have been smoothed using a three point moving average al- 
gorithm. The smoothing technique was applied to each individual /? - curve prior to 
averaging. 

For airfoil, ice shapes and wing geometries the impingement limits on the upper sur- 
face and lower surfaces are given by Sv max and Si maz respectively. For engine inlets the 
impingement limit on the inner is given by S/ mos while that on the outer surface is given by 
So mam - Values for the impingement limits are provided with each averaged impingement 
curve. These values were obtained by comparing the results from the laser reflectance 
method with measurements of the extent of dye on each blotter strip based on visual in- 
spection. The reader is cautioned, however, that the numerical values of the impingement 
limits are approximate because the extent of water impingement cannot be determined 
precisely from the experimental results. The uncertainty in locating the limits of impinge- 
ment is partly due to experimental conditions which varied somewhat between repeated 
runs as discussed above and partly due to the data reduction process. The reflectance 
measurements obtained from the laser reflectometer were normalized by the reflectance 
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value of the bare blotter paper before they were converted into dye mass using the method 
described in section 5.3. Although, the bare blotter paper has nearly uniform reflectance 
properties, small variations do exist. These variations make the prediction of the exact 
locations of the impingement limits difficult. To verify the impingement limits given by the 
laser reflectance method all blotter strips were also inspected visually. The impingement 
limits were obtained by measuring the distance from the highlight to the point where no 
significant amount of dye could be detected by the naked eye. 

The modified inertia parameter, Ko m , d , given in the figures was calculated from the 
following formula: 

*»„, = Ky 

where the inertia parameter, K, was defined in equation (2 — 1) and A/ A s is a function of 
Reynolds number based on MVD and is the ratio of the true range of droplet as projectile 
injected into still air to the range of droplet as projectile following Stoke’s law. The values 
of A/ A * were obtained from the following fourth degree polynomial which is a close fit of 
the curve given in Fig. 6 of Ref. 7. 

— = — 0.022466X 4 + 0.20109X 3 - 0.59067X 2 + 0.36072AT + 0.74544 

A a 


X = log 10 Re mvd , 6 < Re mvd < 1000 


In calculating the inertia parameter, K, the following reference lengths were used: 


1. NLF(l)-0414 airfoil, L = 36 inches 

2. Small and Large Glaze Ice Shapes, L = 2 inches (cylinder diameter) 

3. MS(1)-0317 airfoil, L = 36 inches 

4. NACA 0012 Wing Tip, L = 15/cos(30) = 17.32 inches 

5. Boeing 737-300 Inlet , L = 15 inches ( diameter at compressor face) 


All analysis data presented in this section have been obtained using the experimental 
droplet distributions given in Tables 2.1 and 2.2. Analysis calculations were performed 
with the codes discussed in section 6. 
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7.1 Test Repeatability Results 


Test repeatability for all models tested is demonstrated in Figs. 7.3 through 7.10. 
These figures indicate that in most cases the largest variation in 0 takes place at the point 
of maximum impingement efficiency. The curves presented in Figs. 7.3-7.10 can be used 
to establish error bounds for the experimental data. 

The formula below provides a conservative estimate of test repeatability based on the 
variation in 0max for repeated tests. 


TestRepeatability% 


0.5(0, 


of all run* 


/ max rn i n j rY , um q/ ail r 


unit ^ 


0S 


of all run « 


A percentage of less than 20 % will be regarded as low, and will be considered to be 
indicative of good test repeatability. 


It must be noted that the above equation establishes an upper error bound for the 
experimental data based on a single point on each impingement curve, the point of maxi- 
mum impingement efficiency. In general, all other points on each impingement curve vary 
to a lesser degree from the corresponding points on the average impingement curve. Thus, 
the equation above provides a conservative estimate of test repeatability. 

For the NLF(1)-0414 airfoil test repeatability was approximately ±25% from the av- 
erage impingement curve (Figs. 7.3 and 7.4). The repeatability for the case of a = 8° , 
MVD = 16 fim was ±41% from the average, and indicates that there is a problem with 
the experimental data in this case. 

Test repeatability for the large and small glaze ice shapes was ±31% and ±24% re- 
spectively (see Figs. 7.5 and 7.6). 

Typical repeatability results for the MS(1)~0317 swept airfoil are provided in Fig. 7.7 
and indicate that the individual impingement curves vary by ±17% from the average. 

For the NACA 0012 swept wing tip test repeatability was in the range of ±6% to ±13% 
depending on spanwise location as shown in Fig. 7.8. 


Figures 7,9 and 7.10 present typical test repeatability for the Boeing 737-300 inlet. 
For this geometry test repeatability ranged from ±7% to ±19% depending on MVD and 
circumferential location. 


Impingement limits (“Curve tails”) show good repeatability for all test configurations. 


7.2 Two-Dimensional Impingement Data 


Experimental and analysis water droplet impingement results for the NLF(1)-0414, 
and experimental impingement data for the large and small glaze ice shapes are discussed 
here. The analytical impingement curves for the NLF(1)-0414 airfoil were obtained with 
the computer code discussed in section 6.1. 
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7.2.1 NLF{ 1) -0414 Airfoil Section 

The designation NLF(l)-0414 indicates a natural laminar flow (first) airfoil with a 
design lift coefficient of 0.4 (04) and a thickness of 14 percent chord (Ref. 13). The airfoil 
tested had a 36 inch chord and spanned the height (6 ft) of the tunnel test section. The 
main test variables were angle of attack and spray mean volumetric diameter (MVD). 
Blotter strips were positioned at two spanwise locations (A and B) for each run. Location 
A was 36 inches (midspan) and location B was 32 inches above the tunnel floor. Averaged 
water droplet impingement data for the NLF(l)-0414 at 0 and 8 degrees angle of attack 
and MVDs of 16 and 20pm are summarized in Fig. 7.11. The impingement curves for 
each configuration tested are presented in Figs. 7.12-7.15. 

Referring to Fig. 7.11 the following general trends are observed. For a given angle of 
attack the maximum impingement efficiency, the extent of impingement and the overall 
collection efficiency (area of impingement curve) decrease as MVD is decreased from 20 
pm to 16pm . As the angle of attack is increased, the point of maximum impingement 
moves toward the lower surface and the impingement limits on the lower surface move 
toward the trailing edge. 

The overall agreement between experiment and analysis is good for both MVDs and 
angles of attack. The analytical peak efficiencies, however, are higher at a = 0°, and lower 
at a = 8° than the corresponding experimental values. The main difference between the 
experimental and analytical data is that in the experimental data the maximum impinge- 
ment efficiency for a given MVD increases with angle of attack. The increase in 0 max with 
angle of attack for certain airfoil sections has been observed by other investigators (see 
Ref. 7 pp. 36 and 37). The analytical data for this airfoil, however, show the opposite 
trend. The reason for this discrepancy between test and analysis data is not clear at this 
point. 

7.2.2 Small and Large Glaze Ice Shapes 

Averaged water droplet impingement data for the large and small glaze ice shapes are 
presented in Figs. 7.16 and 7.17 respectively. Both simulated ice shapes had a symmetric 
cross-section and were at zero incidence (i.e., symmetric flow conditions). Thus, symmetric 
0 distributions are expected. The asymmetries observed near the centers of the impinge- 
ment curves for the large and small glaze ice shapes are mainly due to the blotter paper. 
It was very difficult to form the paper to follow precisely the shape of the simulated ice. 
Furthermore, small wrinkles were formed on the blotter strips during their installation 
on the surface of the ice shapes, which modified the reflectance properties of the blotter 
paper. 

The maximum values of impingement efficiency for the large glaze ice shape occur near 
the “horn” regions as shown in Fig. 7.16. These values should be approximately equal 
due to flow symmetry. However, the impingement efficiency at the lower “horn” (see Fig. 

4.2 B, left side looking downstream) has a maximum value of 0.54 while the corresponding 
value on the upper horn is 0.26. Possible reasons for this discrepancy include problems 
with the installation of the blotter strips around the horn regions and flow angularity. 
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7.3 Three-Dimensional Impingement Data 


Experimental water droplet impingement data for the swept MS(1)-0317 airfoil, the 
NACA 0012 swept wing tip and the Boeing 737-300 inlet are presented below. Analytical 
impingement data for the Boeing inlet are also presented and are compared with the 
experimental results. The analysis results for the inlet have been obtained using the 
method described in section 6.2. 

7.3.1 Swept AfS(l) - 317 Airfoil 

This wing had a medium speed (MS) 17% thick airfoil section with a design lift co- 
efficient of 0.3. The wing chord was 36 inches in the streamwise direction and the wing 
span was 6 ft. The wing leading edge sweep angle sweep angle was 30 degrees. This model 
was chosen to investigate the effect of sweep angle on collection efficiency. The model was 
built full span to simulate an infinite swept airfoil. The sweep angle produced a spanwise 
flow thus, the flow was three-dimensional. 

Impingement data were obtained at a = 0 and 8 degrees and for spray clouds with 
MVDs of 16 and 20 pm , Blotter strips were positioned at two spanwise locations A 
and B, 36 and 32 inches above the tunnel floor respectively. The blotter strips were 
in the chordwise (streamwise;) direction. A summary of all experimental water droplet 
impingement data for this swept airfoil is presented in Fig. 7,18. The impingement curves 
for each test case are presented in Figs. 7.19-7.22. 

The experimental results show that the impingement efficiency curves are practically 
the same at the two spanwise locations A and B for all conditions tested. This indicates 
that there was a constant spanwise flow and hence the desired infinite airfoil effect was 
achieved. 

At a = 0 degrees the maximum impingement efficiency is approximately 56% for strip 
A and 59% for strip B for an MVD of 20 pm . The ff max values for the_ 16 pm MVD 
case are 30% for strip A and 39 % for stip B. In ail cases the location of j3 max is on the 
lower surface of the airfoil, near the leading edge. Total water impingement, extent of 
impingement limits and Bmax are less for the smaller MVD size as expected. 

Similar trends are observed at a = 8 degrees with respect to MVD size and strip 
location. For this angle of attack the location of 0 max moves further aft along the lower 
surface and the limits of impingement are significantly greater than for the a = 0 degrees 
case. An interesting feature of the experimental results is that for a given MVD 0 ma x 
increases with angle of attack (see Fig. 7.18). This trend is similar to that observed with 
the NLF(1)-0414 airfoil. 
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7.3.2 NACA 0012 Swept Wing Tip 

Water droplet impingement data for this model are summarized in Fig. 7.23. Detailed 
impingement data are presented in Figs. 7.24-7.27. This was a finite wing with a 30° 
leading edge sweep angle and an aspect ratio of 1.4. The NACA 0012 wing section was 
perpendicular to the leading edge with a chord length of 15 inches. The flow around this 
geometry was highly three-dimensional due to the sweep angle and low aspect ratio of this 
wing. This geometry was tested to investigate the effect of the wing tip vortex on water 
droplet impingement characteristics. 

Tests at a = 0 and 8 degrees were performed for spray clouds with MVDs of 16 and 
20 fim . Blotter strips were positioned at the wing tip, location A, and at two spanwise 
locations, B and C , 42 and 32 inches above the tunnel floor respectively (see Figs. 3.8C 
and 4.4 B). Strip A on the wing tip ran from station B, at the wing leading edge, to about 
50% tip chord (see Fig. 3.8C) The highlight for strip A was located at the leading edge of 
the wing at spanwise location B. There is no upper or lower surface for strip A. Positive 
surface distance for this strip is measured from the leading edge to the trailing edge of the 
wing tip. 

At a = 0° the impingement efficiency curves at spanwise locations B and C are sym- 
metric with respect to the highlight due to symmetric flow conditions. The maximum 
impingement efficiency at these two locations is not significantly affected by MVD or loca- 
tion and it varies from 0.75 to 0.8. The total impingement and the extent of impingement, 
however, is greater for the larger MVD size as expected. The maximum impingement effi- 
ciency for strip A is practically the same as that of strips B and C and is not significantly 
affected by MVD. The total collection efficiency for this strip is also greater for the larger 
MVD. The experimental results show that water impingement at the wing tip was limited 
to the forward 20% of the tip. 

At a -m 8° the point of maximum impingement and the impingement limits at locations 
B and C shift slightly toward the lower surface. The maximum impingement efficiency, 
total water droplet impingement and extent of impingement limits are greater for the larger 
MVD. Similar trends with MVD are observed for strip A. The extent of impingement for 
strip A is considerably reduced at the higher angle of attack due to the presence of the tip 
vortex. As angle of attack is increased the vortex strength and vortex core diameter are 
increased. The larger vortex causes a larger disturbance to the flow locally and a reduction 
in water droplet impingement. 

Referring to Fig. 7.23 (summary of results) certain trends in the experimental data 
require further discussion. For the 20 pm MVD case, the maximum impingement efficiency 
at locations B and C are very similar for both angles of attack. This is not the case for 
strip A for which impingement efficiency is higher at a = 8° . It is not clear why 0 m ax 
for strip A should increase with Zof attack. For MVD = 16 fim , the results show that 
0max a = 8° at all spanwise locations is lower than the corresponding fimax values 
at a = 0° . Furthermore, 0 max in this case is a function of location, being maximum at 
location A and minimum at location C. The reasons for the observed differences in /? max 
for the two MVD sizes are not clear. It is possible that these differences in trends axe 
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due to experimental error in achieving the desired experimental conditions. However, test 
repeatability for this test geometry at a = 8° was very good. 

A few comments regarding the impingement limits are also in order. The experimental 
impingement curves for 0 = 0°, indicate that the extent of impingement is not significantly 
affected by MVD. At a — 8° several of the impingement curves fall off abruptly. These 
discrepancies are attributed to the reflectance properties of the blotter paper which were 
modified by the adhesive tape used to mount the strips on the wing. The correct values for 
the impingement limits are the numerical values given below each curve. These values were 
obtained by careful examination of the raw laser reflectance data and by visual inspection 
of the blotter strips. 


7.3.3 Boeing 737-300 Engine Inlet 

This inlet had a non-axisymmetric three-dimensional geometry which, when placed 
at a high angle of attack, produced a complex flow pattern. The flowfield was compli- 
cated further by variations in mass flow through the inlet used to simulate different inlet 
capture area ratios. The resulting flow configuration was a challenging test case for the 
experimental and analytical methods. 

For this geometry, the test variables included a , MVD and inlet mass flow. The high 
mass flow (23 lbm/sec) corresponds to an inlet capture area ratio of 1. The low mass flow 
(17 lbm/sec) corresponds to a capture area ratio of 0.75. 

A total of 8 blotter strips were positioned around the inlet for each test. The strips 
were set at 45 degrees intervals corresponding to 0 = 0, 45, 90, 135, 180, 225, 270 and 325 
degrees (see Fig. 3.7B). Two to five runs were made for each combination of a , MVD 
and mass flow setting. Results from several of these runs were averaged to produce the 
impingement curves shown in Figs. 7.28-7.30. Analytical data are also presented in these 
figures. Figure 7.28 provides a summary of all experimental results obtained for the Boeing 
737-300. Detailed results for the 0 = 0° and a = 15° test cases are presented in Fig. 7.29 
and Fig. 7.30 respectively. 

The engine inlet model was mounted in the tunnel with the top (0 = 0° ) and bottom 
(0 = 180° ) of the inlet in a horizontal plane. Due to the symmetry of the flow with respect 
to a horizontal plane through locations A ( 0 ~ 0° ) and E ( 0 = 180° ) the impingement 
curves for strips B and H should be the same. This is also true for strips C and G and 
strips D and F. Thus, all B and H strips for a given test condition were combined into a 
single averaged impingement curve. The same was done with all C and G and all D and 
F strips. 

Water droplet impingement data for this inlet were presented in Reference 10. These 
data were obtained in the September 1985 tests. Discrepancies between the analytical and 
the 1985 test results for the case of a = 15 degrees were observed and these differences 
(particularly for 0 = 0 degrees) were attributed to non-uniform dye distribution in the 
droplet cloud. To clarify some of these discrepancies some of the test points were repeated 
during the current tests. Extra attention was paid to the cloud uniformity in the current 
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tests. The result being that the new experimental data are in much better agreement with 
the analytical data. 

At a = 0 degrees the following trends are observed from the experimental data pre- 
sented in Fig. 7.29: 


a . Reducing the inlet mass flow from 23 to 17 Ibm/s causes the location of /3 moz 
to shift toward the inner cowl (positive s). 

b . For the same mass flow and cloud MVD, the location of 0 max moves progres- 
sively toward the inner cowl as 0 varies from 0 to 180 degrees. The maximum 
impingement limits are observed at 0 = 135° and 0 = 180° . 

c . For the low mass flow and cloud MVD the area under the impingement curves 
is reduced indicating lower total water impingement. 

d , The analytical impingement curves are wider than the experimental curves. 

The correlation between experiment and theory is good in most cases. The 
analysis data indicate that the maximum value of j3 max occurs at 0 — 90° 
in all cases presented and that the maximum impingement limits are at 
0 = 135° . 

Referring to Fig. 7.30 the following trends are observed for the case a — 15 degrees. 

a . For a given circumferential location (constant 0 )_ reducing the inlet mass 
flow from 23 to 17 Ibm/s causes the location of /? maz to shift toward the 
inner cowl (positive s) 

b . For the same mass flow and cloud MVD, the location of /f maz moves pro- 
gressively from the inner cowl at 0 = 0° to the outer cowl at 0 = 180° . The 
maximum impingement limits are observed at 0 =135° 

c . For the low mass flow and cloud MVD the area under the impingement curves 
is reduced indicating lower total water impingement. 

d . The analytical impingement curves are in all cases wider than the experi- 
mental curves. The correlation between experiment and theory is good in 
most cases. The experimental values of 0 m ax are in general higher than the 
corresponding analytical values. The analysis results indicate that the extent 
of the impingement limits is maximum at 0 = 135° and the maximum 0 max 
occurs 0 = 90° . 

More details regarding impingement trends for the Boeing 737-300 inlet can be found 
in Ref. 10. 
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TABLE 7.1 

REFERENCE COLLECTOR DYE MASS FOR AIRFOILS, 
ICE SHAPES AND WING TIP 


Test Condition 

Collector Dye Mass in fig /cm 2 

a “ 0, 8° , MVD = 20 fim , Locations A and B 

2.223 

a — 0,8°, MVD = 16 fim , Locations A and B 

2,608 


TABLE 7.2 

REFERENCE COLLECTOR DYE MASS FOR ENGINE INLET 


Test Condition 

Collector Dye Mass in figjcrn* 

$ = 0° 
Strip A 

$ m 45° 
Strip B 

$ = 90° 
Strip C 

9 = 135° 
Strip D 

B = 180° 
Strip E 

8 = 225° 
Strip F 

6 = 270° 
Strip G 

9 = 315° 
Strip H 

a = 0°, MVD= 20 fim 

2.1738 

2.0456 

2,0968 

2.3740 

2.5218 

2.3638 

2.2723 

1.8466 

a = 0°j MVD— 16 fim 

2.4418 

2.2987 

2.4401 

2.5842 

2.6127 

2*4975 

2.5977 

2.1527 

a = 15°, MVD= 20 ftm 

1.5618 

1.5514 

1.8894 

2.0640 

2.2003 

2.0755 

1.7509 

1.6261 

a = 15° , MVD= 16 fim 

2.2866 

2.3565 

2.5111 

2.6744 

2.6926 

2.5766 

2.6116 

2.3594 
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FIGURE 7.1 

HIGHLIGHT MARKS AND SURFACE DISTANCE CONVENTION 
FOR WINGS AND ICE SHAPES 




TEST RUN ID: 220N900A, 221N900A, 222N900A, 
223N900A, 224N900A, 225N900A 



(Positive Surface Distance, S, on Lower Surface) 

(A) MVD = 20 MICRONS, STRIP A 
FIGURE 7.3 

TEST REPEATABILITY FOR NLF(l)-0414 AIRFOIL AT « = 0° 

(PAGE 1 OF 2). 
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LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 226N600A, 227N600A, 228N60QA, 

229N600A, 230N600A 



(Positive Surface Distance, S, on Lower Surface) 

(B) MVD = 16 MICRONS, STRIP A 
FIGURE 7.3 

TEST REPEATABILITY FOR NLF(1)-0414 AIRFOIL AT a = 0° 

(PAGE 2 OF 2). 
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LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 215N908A, 216N908A, 217N908A, 

218N908A, 219N908A 



(Positive Surface Distance, S, on Lower Surface) 
(A) MVD = 20 MICRONS, STRIP A 


FIGURE 7.4 

TEST REPEATABILITY FOR NLF(1)~0414 AIRFOIL AT « = 8° 

(PAGE 1 OF 2). 
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LOCAL IMPINGEMENT EEEIOIE2ST 


TEST RUN ID: 232N608A, 233N608A, 234N608A, 

235N608A, 236N608A 



(Positive Surface Distance, S, on Lower Surface) 

(B) MVD = 16 MICRONS, STRIP A 
FIGURE 7.4 

TEST REPEATABILITY FOR NLF(1)-0414 AIRFOIL AT a = 8° 

(PAGE 2 OF 2). 
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1.00 


TEST RUN ID: 121I900A, 122I900A, 123I900A, 
124I900A, I25I900A, 1261900 A 


0 . 00 - 

* 

Do.80 

z 

« 

W 

0 0 . 70 - : 

w 

h 



SURFACE DISTANCE FROM HIGHLIGHT - mm 

i 


(Positive Surface Distance, S, on Lower Surface) 

MVD = 20 MICRONS, STRIP A 
FIGURE 7.6 

TEST REPEATABILITY FOR SMALL GLAZE ICE SHAPE AT a = 0° 
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OOl 



(Positive Surface Distance, S, on Lower Surface) 

MVD = 20 MICRONS, STRIP A 
FIGURE 7.7 

TEST REPEATABILITY FOR MS(1)-0317 AIRFOIL AT a = 8° 



180 
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TEST RUN ID: 127W900A, 128W900A, 129W900A 

130W900A, 131W900A 


SURFACE DISTANCE FROM HIGHLIGHT — xxxzxt 
(Positive Surface Distance, S, on Lower Surface) 

(A) MVD = 20 MICRONS, STRIP A 
FIGURE 7.8 

TEST REPEATABILITY FOR NACA 0012 SWEPT WING TIP AT « 

(PAGE 1 OF 3). 




LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 127W900B, 128W900B, 129W900B, 

130W900B, I3IW900B 



(Positive Surface Distance, S, on Lower Surface) 

(B) MVD = 20 MICRONS, STRIP B 
FIGURE 7.8 

TEST REPEATABILITY FOR NACA 0012 SWEPT WING TIP AT a = 0° 

(PAGE 2 OF 3). 
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LO CA1L IMPINGEMENT EFFICIEN 


TEST RUN ID; 127W900C, 128W900C, 129W900C, 

130W900C, 131W900C 



(C) MVD = 20 MICRONS, STRIP C 


FIGURE 7.8 

TEST REPEATABILITY FOR NACA 0012 SWEPT WING TIP AT a = 0° 

(PAGE 3 OF 3). 
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TEST RUN ID: 17439I5DF, 17539I5DF, 1763915DF 



SURFACE DISTANCE FROM HIGHLIGHT — mm 


(Positive Surface Distance, S, on Inner Surface) 

(B) MVD = 20 MICRONS, MASS FLOW = 22.96 LBM/S, 6 = 135° , 225° 

FIGURE 7.9 

TEST REPEATABILITY FOR BOEING 737-300 AIRFOIL AT a = 15° 

(PAGE 2 OF 3). 
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LOCAL HVdCJPUST GEMENT EFFICIENCY 


TEST RUN ID: 1743915E, 1753915E, 1763915E 



(Positive Surface Distance, S, on Inner Surface) 

(G) MVD = 20 MICRONS, MASS FLOW = 22,96 LBM/S, 0 = 180° 

FIGURE 7.9 

TEST REPEATABILITY FOR BOEING 737-300 INLET AT a = 15° 

(PAGE 3 OF 3). 
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LOCAL IMPINGEMENT EFFICIEN 


TEST RUN ID; 1683615A, 1693615A, 1703615A, 
1713615A, 1723615A, 1733615A 



(Positive Surface Distance, S, on Inner Surface) 

(A) MVD = 16 MICRONS, MASS FLOW = 17.2 LBM/S, 0 = 0 ° 

FIGURE 7.10 

TEST REPEATABILITY FOR BOEING 737-300 INLET AT a = 15° 

(PAGE 1 OF 5). 


78 



LOCAL IMPINGEMENT EFFICIENCY 


I.GCh 


TEST RUN ID: 1683615BH, 1693615BH, 1703615BH 


0.9CH 


O.SOI 



(Positive Surface Distance, S, on Inner Surface) 


(B) MVD = 16 MICRONS, MASS FLOW = 17.2 LBM/S, 6 = 45° , 315° 

FIGURE 7.10 

TEST REPEATABILITY FOR BOEING 737-300 INLET AT a = 15° 

(PAGE 2 OF 5). 
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LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 1683615CG, I693615CG, 1703615CG 



(Positive Surface Distance, S, on Inner Surface) 

(C) MVD = 16 MICRONS, MASS FLOW = 17.2 LBM/S, 9 = 90° , 270° 

FIGURE 7.10 

TEST REPEATABILITY FOR BOEING 737-300 INLET AT a = 15° 

(PAGE 3 OF 5). 
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TEST RUN ID: 1683615DF, 1693615DF, I703615DF 


SURFACE DISTANCE FROM HIGHLIGHT — mm 
(Positive Surface Distance, S, on Inner Surface) 

(D) MVD = 16 MICRONS, MASS FLOW = 17.2 LBM/S, 0 = 135 

FIGURE 7.10 

TEST REPEATABILITY FOR BOEING 737-300 INLET AT a = 

(PAGE 4 OF 5). 




LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 1683615E, 1693615E, 17Q3615E, 
1713615E, 1723615E, 1733615E 



(Positive Surface Distance, S, on Inner Surface) 

(E) MVD = 16 MICRONS, MASS FLOW = 17.2 LBM/S, 0 = 180° 

FIGURE 7.10 

TEST REPEATABILITY FOR BOEING 737-300 INLET AT a = 15° 

(PAGE 5 OF 5). 
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LOCAL IMPINGEMENT EFFICIENCY { $ ) - DIMENSIONLESS 











LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 220N900A, 221N900A, 222N900A, 
223N900A, 224N900A, 225N900A 



(Positive Surface Distance, S, on Lower Surface) 
ffimaz = 0,373, Su max = —48 mm, Si, max = 28 mm (Experiment) 

(A) STRIP A, MVD = 20 MICRONS 

FIGURE 7.12 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR NLF(1)-0414 AIRFOIL AT a = 0° (^=165 mph, K 0 d = 0.036) 

(PAGE 1 OF 2). 
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TEST RUN ID: 220N900B, 221N900B, 222N900B, 
223N900B, 224N900B, 225N900B 



(Positive Surface Distance, S, on Lower Surface) 
fima* = 0.372, Su max = — 48mm, Si, max = 28mm (Experiment) 

(A) STRIP B, MVD = 20 MICRONS 

FIGURE 7.12 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR NLF(l)~0414 AIRFOIL AT a = 0° (^=165 mph, K 0 „ = 0.036) 

(PAGE 2 OF 2). 
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TEST RUN ID: 226N600A, 227N600A, 228N600A, 

229N600A, 230N600A 



(Positive Surface Distance, S, on Lower Surface) 

0max — 0.344, Su maK — — 37mm, = 20mm (Experiment) 

(A) STRIP A, MVD = 16 MICRONS 

FIGURE 7.13 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR NLF(l)-G414 AIRFOIL AT a = 0° (^=165 mph, K 0 d = 0.025) 

(PAGE 1 OF 2). 
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TEST RUN ID: 226N600B, 227N600B, 228N600B, 

229N600B, 230N600B 



(Positive Surface Distance, S, on Lower Surface) 

0max = 0.285, $v max = —37 mm, Si, mam = 20 mm (Experiment) 

(A) STRIP B, MVD = 16 MICRONS 

FIGURE 7.13 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR NLF(1)-0414 AIRFOIL AT a = 0° (^=165 mph, K 0mvd = 0.025) 

(PAGE 2 OF 2). 
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LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 215N908A, 216N908A, 2I7N908A, 

218N908A, 219N908A 



(Positive Surface Distance, S, on Lower Surface) 

Pmax = 0.562, Su max = -10mm, Sjr, moi = 230mm (Experiment) 

(A) STRIP A, MVD = 20 MICRONS 

FIGURE 7.14 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR NLF(1)-0414 AIRFOIL AT a = 8° (Voo=165 mph, K 0m , d = 0.036) 

(PAGE 1 OF 2). 
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LO CAL IMPINGEMENT EF’FIOIEISrCY 


TEST RUN ID: 215N908B, 216N908B, 217N908B, 

218N908B, 219N908B 



(Positive Surface Distance, S, on Lower Surface) 

/ 3 mo * = 0.566, S Urnai = -10 mm, S Lmax = 230mm (Experiment) 

(A) STRIP B, MVD = 20 MICRONS 

FIGURE 7.14 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR NLF(l)-0414 AIRFOIL AT a = 8° (Voo=165 mph, K 0<n , d = 0.036) 

(PAGE 2 OF 2). 
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L.OCAJL. IMPINGEMENT EFFICIENCY 


TEST RUN ID: 232N608A, 233N608A, 234N608A, 

235N608A, 236N608A 



(Positive Surface Distance, S, on Lower Surface) 

0max = 0.382, Sv max = —5 mm, Si, m « = 127mm (Experiment) 

(A) STRIP A, MVD = 16 MICRONS 

FIGURE 7.15 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR NLF(1)-0414 AIRFOIL AT a = 8° (Voo=165 mph, K 0m d = 0.025) 

(PAGE 1 OF 2). 
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local impingement efficiency 


TEST RUN ID: 232N608B, 233N608B, 234N608B, 

235N608B, 236N608B 



(Positive Surface Distance, S, on Lower Surface) 

Pmax = 0,368, = —5mm, Si, max = 127mm (Experiment) 

(A) STRIP B, MVD = 16 MICRONS 

FIGURE 7.15 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR NLF(l)-0414 AIRFOIL AT a = 8° (1^=165 mph, K 0mvd = 0.025) 

(PAGE 2 OF 2). 
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LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 116I900AB, 117I900AB, 11 81900 AB, 

1 191900 AB, 1 201900 AB 



(Positive Surface Distance, S, on Lower Surface) 

0max = 0.540, Su nma = -53 mm, S Lmetx = 53mm (Experiment) 

(A) STRIPS A,B, MVD = 20 MICRONS 

FIGURE 7.16 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR LARGE GLAZE ICE SHAPE, MOUNTED ON 2 INCH CYLINDER, 
AT a = 0° (Voo=165 mph, K 0m , d = 0.648). 
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LOCAL IMPIN CEMENT EFFICIENCY 


TEST RUN ID: 121I900A, 122I900A, 123I900A, 

124I900A, 125I900A 



(Positive Surface Distance, S, on Lower Surface) 

0m*x = 0.460, S Umaz = "37mm, S Lrnax = 37mm (Experiment) 

(A) STRIPS A,B, MVD = 20 MICRONS 

FIGURE 7.17 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR SMALL GLAZE ICE SHAPE, MOUNTED ON 2 INCH CYLINDER, 
AT et = 0° (V 00 =: 165 mph, K 0mvd = 0.648). 
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LOCAL IMPINGEMENT EFFICIENCY ( f} ) - DIMENSIONLESS 
















LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 149M900A, 150M900A, 151M900A, 

152M900A, 153M900A 



(Positive Surface Distance, S, on Lower Surface) 

= 0.558, S Ummm = -63mm, S Lmax = 65 mm 

(A) STRIP A, MVD = 20 MICRONS 

FIGURE 7.19 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR MS(1)-0317 AIRFOIL AT « = 0° (¥^=165 mph, K 0mvd = 0.036) 

(PAGE 1 OF 2). 
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LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 149M900B, 150M900B, 151M900B, 

152M900B, 153M900B 



(Positive Surface Distance, S, on Lower Surface) 

0max = 0,593, ~ 63mm, Si Amax — 65 mm 

(B) STRIP B, MVD = 20 MICRONS 

FIGURE 7.19 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR MS(l)~0317 AIRFOIL AT a — 0° (^=165 mph, K 0mvd = 0.036) 

(PAGE 2 OF 2). 
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LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 154M600A, 155M600A, 156M600A, 

157M600A 



(Positive Surface Distance, S, on Lower Surface) 
fimaz = 0.302, Sv max = —40 mm, Si, max = 30mm 

(A) STRIP A, MVD = 16 MICRONS 

FIGURE 7.20 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR MS(1)-0317 AIRFOIL AT a = 0° (^=165 mph, K 0mtd = 0.025) 

(PAGE 1 OF 2). 
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LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 154M600B, 155M600B, 156M600B, 

157M600B 



(Positive Surface Distance, S, on Lower Surface) 
fimaz = 0.395, Su max = -40mm, Sj, mM = 29mm 

(B) STRIP B, MVD = 16 MICRONS 

FIGURE 7.20 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR MS(1)-0317 AIRFOIL AT a = 0° (^=165 mph, K 0mad = 0.025) 

(PAGE 2 OF 2). 
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LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 163M908A, 164M908A, 165M908A, 

166M908A, 167M908A 



(Positive Surface Distance, S, on Lower Surface) 

0max = 0.595, Su m9x = —12 mm, S Lmai = 155mm 

(A) STRIP A, MVD = 20 MICRONS 

FIGURE 7.21 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR MS(1)-0317 AIRFOIL AT a =.8° {^=165 mph, K 0mvd = 0.036) 

(PAGE 1 OF 2). 
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LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 163M908B, 164M908B, 165M908B, 

166M908B, I67M908B 



(Positive Surface Distance, S, on Lower Surface) 

Pmax = 0.597, Su m „ — -10mm, Sx, mai = 155mm 

(B) STRIP B, MVD = 20 MICRONS 

FIGURE 7.21 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR MS(1)-0317 AIRFOIL AT a = 8° (^=165 mph, K 0mvi = 0.036) 

(PAGE 2 OF 2). 
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LOCAL IMPINGEMENT EFPICIEN 


TEST RUN ID: 158M608A, 159M608A, 160M608A, 

161M608A, 162M608A 



(Positive Surface Distance, S, on Lower Surface) 

Ima* = 0.475, S Umax = -4mm, S Lmax - 106mm 

(A) STRIP A, MVD = 16 MICRONS 

FIGURE 7.22 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR MS(1)-0317 AIRFOIL AT « = 8° (1^=165 mph, K 0m , d = 0.025) 

(PAGE- 1 OF 2). 
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LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 158M608B, I59M608B, 160M608B, 

161M608B, 162M608B 



(Positive Surface Distance, S, on Lower Surface) 
fimax — 0,472, Su max = —4mm, Si, maz = 110mm 

(B) STRIP B, MVD = 16 MICRONS 

FIGURE 7.22 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR MS(1)-0317 AIRFOIL AT a = 8° (V 0O =165 mph, K 0m , d = 0.025) 

(PAGE 2 OF 2). 
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LOCAL IMPINGEMENT EFFICIENCY ( £ } - DIMENSIONLESS 







LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 127W900A, 128W900A, 129W900A, 

130W900A, 131W900A 



(Positive Surface Distance, S, on Lower Surface) 

0max = 0.774, Stotai =: 76 mm 

(A) STRIP A, MVD = 20 MICRONS 

FIGURE 7.24 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR NACA 0012 SWEPT WING TIP AT a = 0° (^=165 mph, K 0mvd = 0.075) 

(PAGE 1 OF 3). 
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TEST RUN ID: 127W9G0B, 128W9Q0B, 129W900B, 

130W900B, 131W900B 



(Positive Surface Distance, S, on Lower Surface) 
ftmax — 0.750, Su max — —40 mm, St max — 40 mm 

(B) STRIP B, MVD = 20 MICRONS 

FIGURE 7.24 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR NACA 0012 SWEPT WING TIP AT« = 0° (^<*=165 mph, K 0 . = 0.075) 

(PAGE 2 OF 3). 
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LOCAL impingement efficiency 


TEST RUN ID: 127W900C, 128W900C, I29W900C, 

130W900C, 131W9G0C 



(Positive Surface Distance, S, on Lower Surface) 

0max = 0.802, Su max = — 40mm, = 40mm 

(C) STRIP C, MVD = 20 MICRONS 

FIGURE 7.24 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR NACA 0012 SWEPT WING TIP AT a = 0° (^=165 mph, K 0m , d = 0.075) 

(PAGE 3 OF 3). 
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LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 132W600A, 133W600A, 134W600A 

135W600A, 136W600A 


Q TEST DATA 



•120 -90 -60 -30 0 30 

SURFACE DISTANCE FROM HIGHLIGHT 


mm 


(Positive Surface Distance, S, on Lower Surface) 

0max = 0.798, Stotai = 64mm 

(A) STRIP A, MVD = 16 MICRONS 

FIGURE 7.25 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR NACA 0012 SWEPT WING TIP AT « = 0° (Voo=165 mph, K 0m , d = < 

(PAGE 1 OF 3). 
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LOCAL IMF* IN GEMENT EFFICIENCY 


TEST RUN ID: 132W600B, 133W600B, 134W600B, 

135W600B, 136W600B 



(Positive Surface Distance, S, on Lower Surface) 
fimax = 0,773, Sv mai = -30mm, S Lmax = 30mm 

(B) STRIP B, MVD = 16 MICRONS 

FIGURE 7.25 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR NACA 0012 SWEPT WING TIP AT a = 0° (^=165 mph, K 0mvd = 0.052) 

(PAGE 2 OF 3). 
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TEST RUN ID: 132W600C, 133W600C, 134W600C, 

135W600C, 136W600C 



(Positive Surface Distance, S, on Lower Surface) 

0max — 0.755, S Vmax = -30 mm, S Lmax = 30mm 

(C) STRIP C, MVD = 16 MICRONS 

FIGURE 7.25 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR NACA 0012 SWEPT WING TIP AT a = 0* (^=165 mph, K 0m , d = 0.052) 

(PAGE 3 OF 3). 
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LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 143W908A, 144W908A, 145W908A, 

146W908A, 147W908A 



(Positive Surface Distance, S, on Lower Surface) 

~~~ 0.889, S total “ 50mm 

(A) STRIP A, MVD = 20 MICRONS 

FIGURE 7.26 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR NACA 0012 SWEPT WING TIP AT a = 8° (Voo=165 mph, K 0mvd = 0.075) 

(PAGE 1 OF 3). 
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TEST RUN ID: 143W908B, 144W908B, 145W9G8B 

146W908B, 147W908B 


TEST DATA 


SURFACE DISTANCE FROM HIGHLIGHT — mm 

(Positive Surface Distance, S, on Lower Surface) 
flmax = 0.790, S{j mas = -12mm, S Ltnax = 100mm 

(B) STRIP B, MVD = 20 MICRONS 

FIGURE 7.26 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR NACA 0012 SWEPT WING TIP AT a = 8° (F 00 =165 mph, K 0mvd = ( 

(PAGE 2 OF 3). 




LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 143W908C, 144W908C, 145W908C, 

146W908C, 147W908C 



(Positive Surface Distance, S, on Lower Surface) 

$max = 0.770, S Umai = -13mm, S Lmair = 101mm 

(C) STRIP C, MVD = 20 MICRONS 

FIGURE 7.26 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR NACA 0012 SWEPT WING TIP AT a = S° (Voo=165 mph, K 0mvd = 0.075) 

(PAGE 3 OF 3). 
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LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 137W608A, 138W608A, 139W608A, 

140W608A, 141W608A 



(Positive Surface Distance, S, on Lower Surface) 

$max = 0,721, Stotai = 49mm 

(A) STRIP A, MVD = 16 MICRONS 

FIGURE 7.27 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR NACA 0012 SWEPT WING TIP AT a = 8° (^*>=165 mph, K 0 d = 0.052) 

(PAGE 1 OF 3). 
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LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 137W608B, 138W608B, 139W608B, 

140W608B, 141W608B 



(Positive Surface Distance, S, on Lower Surface) 

0max = 0.650, Sy max = —9 mm, Sx m<ut = 67 mm 

(B) STRIP B, MVD = 16 MICRONS 

FIGURE 7.27 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR NACA 0012 SWEPT WING TIP AT a = 8° (^*>=165 mph, K 0 d = 0.052) 

(PAGE 2 OF 3). 
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LOCAL IMPINGEMENT EF^ICIEISTOY 


TEST RUN ID: 137W608C, 138W608C, 139W608C, 

140W608C, 141W608C 



(Positive Surface Distance, S, on Lower Surface) 

Pma* = 0.557, S Um „ = -9mm, S Lmma = 67 mm 

(C) STRIP C, MVD = 16 MICRONS 

FIGURE 7.27 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR NACA 0012 SWEPT WING TIP AT a = 8° (^=165 mph, K 0mmd = 0.052) 

(PAGE 3 OF 3). 
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LOCAL IMPINGEMENT EFFICIENCY ( fi ) - DIMENSIONLESS 


9 = 0 ° 


5 = 45° 


5*60° 


$ ta 135 * 


9 = 180 * 



W=23 

MVD=20 

pr=0° 

W=17 

MVD=16 

oeO’ 

W«» 

MVD-20 

o»lS* 

W=17 

MVD»16 

owls’ 


-120 0 120 -120 0 120 -120 0 120 -120 0 120 -120 0 120 


SURFACE DISTANCE FROM HIGHLIGHT - MM 

NOTE : W - LBM / SEC, MVD - pm 

NEGATIVE SURFACE DISTANCE - OUTER SURFACE 
POSITIVE SURFACE DISTANCE - INNER SURFACE 

FIGURE 7.28 



AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR BOEING 737-300 INLET. SUMMARY OF RESULTS. 




LOCAL XMPXCTGEKfENT EFFICIENCY 


TEST RUN ID: 1773900A, 1783900A 



(Positive Surface Distance, S, on Inner Surface) 

0 ma x = 0.440, So max — -60mm, S Jmaz = 20mm ( Experiment ) 

(Al) STRIP A, 8 = 0° , MVD = 20 MICRONS, MASS FLOW = 23 LBM/SEC 

FIGURE 7.29 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR BOEING 737-300 INLET AT a = 0° (^=173 mph, K 0mvd = 0.089) 

(PAGE I OF 10). 
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LOCAL IMPINGEMENT EF'F'ICIEINrC'V 


TEST RUN ID: 1773900BH, 1783900BH 



(Positive Surface Distance, S, on Inner Surface) 

0max — 0.680, $o maat = — 55mm, Sj max = 23 mm ( Experiment ) 

(A2) STRIPS B,H, 0 = 45° , MVD = 20 MICRONS, MASS FLOW = 23 LBM/SEC 

FIGURE 7.29 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR BOEING 737-300 INLET AT a = 0° (^=173 mph, K 0mvd = 0.089) 

(PAGE 2 OF 10). 
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TEST RUN ID: 1773900CG, 1783900CG 



(Positive Surface Distance, S, on Inner Surface) 

~ 0.540, So max = — 43mm, = 28mm ( Experiment ) 

(A3) STRIPS C,G, $ = 90° , MVD = 20 MICRONS, MASS FLOW = 23 LBM/SEC 

FIGURE 7.29 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR BOEING 737-300 INLET AT a = 0° (^=173 mph, K 0 , = 0.089) 

(PAGE 3 OF 10). 
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LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID; 1773900DF, 1783900DF 



(Positive Surface Distance, S, on Inner Surface) 

Pmax = 0.520, So max = — 46mm, Sj max — 34mm ( Experiment ) 

(A4) STRIPS D,F, 6 = 135° , MVD = 20 MICRONS, MASS FLOW = 23 LBM/SEC 

FIGURE 7.29 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR BOEING 737-300 INLET AT a = 0° (1^=173 mph, K 0 vd = 0.089) 

(PAGE 4 OF 10). 
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TEST RUN ID: 1773900E, 1783900E 



(Positive Surface Distance, S, on Inner Surface) 
ft max = 0.560, So max = — 38mm, S/ max = 40mm ( Experiment ) 

(A5) STRIP E, 9 = 180° , MVD = 20 MICRONS, MASS FLOW = 23 LBM/SEC 

FIGURE 7.29 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR BOEING 737-300 INLET AT a = 0° (^=173 mph, K 0mvd = 0.089) 

(PAGE 5 OF 10). 
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LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 1793600A, 18 03 600 A 



(Positive Surface Distance, S, on Inner Surface) 

Pmax = 0,550, So maz = -28mm, Sj maK = 23mm ( Experiment ) 

(Bl) STRIP A, 0 = 0° , MVD = 16 MICRONS, MASS FLOW = 17 LBM/SEC 

FIGURE 7,29 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR BOEING 737-300 INLET AT a = 0° (^=173 mph, K 0mvd = 0.062) 

(PAGE 6 OF 10). 
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TEST RUN ID: 1793600BH, 1803600BH 



(Positive Surface Distance, S, on Inner Surface) 

0max ^ 0.620, S 0max - -21mm, S Jmmm = 23mm ( Experiment ) 

(B2) STRIPS B,H, $ = 45° , MVD = 16 MICRONS, MASS FLOW = 17 LBM/SEC 

FIGURE 7.29 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR BOEING 737-300 INLET AT a = 0° (^=173 mph, K 0mnd = 0.062) 

(PAGE 7 OF 10). 
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LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 1793600CG, 1803600CG 



(Positive Surface Distance, S, on Inner Surface) 

0max = 0.570, So max = -16mm, Sj max = 29mm ( Experiment ) 

(B3) STRIPS C,G, 6 = 90° , MVD = 16 MICRONS, MASS FLOW = 17 LBM/SEC 

FIGURE 7.29 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR BOEING 737-300 INLET AT a = 0° (Voo=173 mph, K 0m , d = 0.062) 

(PAGE 8 OF 10). 
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LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 1793600DF, 18036Q0DF 



(Positive Surface Distance, S, on Inner Surface) 

Pmax = 0.570, S 0maz = -20mm, S tmaM = 36mm ( Experiment ) 

(B4) STRIPS D,F, 0 = 135° , MVD = 16 MICRONS, MASS FLOW = 17 LBM/SEC 

FIGURE 7.29 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR BOEING 737-300 INLET AT a = 0° (Voo=173 mph, K 0m , d = 0.062) 

(PAGE 9 OF 10). 



LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 1793600E, 1803600E 



(Positive Surface Distance, S, on Inner Surface) 

Pmax = 0.540, So max = —15mm, Sj max = 40 mm ( Experiment ) 

(B5) STRIP E, 8 = 180° , MVD = 16 MICRONS, MASS FLOW = 17 LBM/SEC 

FIGURE 7.29 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR BOEING 737-300 INLET AT a = 0° (^=173 mph, K 0 , d = 0.062) 

(PAGE 10 OF 10). 
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LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 1743915 A, 1753915 A, 1763915A 



(Positive Surface Distance, S, on Inner Surface) 
fimax = 0,710, So mat = — 22mm , Si max = 43 ram ( Experiment ) 

(Al) STRIP A, 8 = 0° , MVD = 20 MICRONS, MASS FLOW = 23 LBM/SEC 

FIGURE 7.30 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR BOEING 737-300 INLET AT a = 15 c (Voo-171 mph, K 0mvd = 0.089) 

(PAGE 1 OF 10). 
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LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 1743915BH, 1753915BH, 1763915BH 



(Positive Surface Distance, S, on Inner Surface) 

Pmax = 0.740, S 0max = ~25mm, Sr max = 39 mm ( Experiment ) 

(A2) STRIPS B,H, $ = 45° , MVD = 20 MICRONS, MASS FLOW = 23 LBM/SEC 

FIGURE 7.30 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR BOEING 737-300 INLET AT a = 15° (Voo=171 mph, K 0mvd = 0.089) 

(PAGE 2 OF 10). 
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LOOAJL IMPINGEMENT EFFICIEN CY 


TEST RUN ID: 1743915CG, 1753915CG, 1763915CG 



(Positive Surface Distance, S, on Inner Surface) 
flmax ~ 0.650, S 0maz = -39mm, S [max = 28mm ( Experiment ) 

(A3) STRIPS C,G, 0 = 90° , MVD = 20 MICRONS, MASS FLOW = 23 LBM/SEC 

FIGURE 7.30 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR BOEING 737-300 INLET AT a = 15° (V TO =171 mph, K 0mvd = 0.089) 

(PAGE 3 OF 10). 
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LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 1743915DF, 1753915DF, 1763915DF 



(Positive Surface Distance, S, on Inner Surface) 

0max = 0.600, So max = — 190mm, Sj mmm = 15mm ( Experiment ) 

(A4) STRIPS D,F, $ = 135° , MVD = 20 MICRONS, MASS FLOW = 23 LBM/SEC 

FIGURE 7.30 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR BOEING 737-300 INLET AT a = 15° (^=171 mph, K 0 d = 0.089) 

(PAGE 4 OF 10). 
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LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 1743915E, 17539I5E, 1763915E 



(Positive Surface Distance, S, on Inner Surface) 
j 9maz = 0.510, S 0max = -187mm, S Jmai = 13mm ( Experiment ) 

(A5) STRIP E, 9 = 180° , MVD = 20 MICRONS, MASS FLOW = 23 LBM/SEC 

FIGURE 7.30 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR BOEING 737-300 INLET AT a = 15° (Voo=171 mph, K 0mvd = 0.089) 

(PAGE 5 OF 10). 


131 



LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 1683615A, 1693615A, 1703615A, 1713615A, 1723615A, 17336I5A 



(Positive Surface Distance, S, on Inner Surface) 

Pmax = 0.700, So max = — 10mm, Sj max = 48mm ( Experiment ) 

(Bl) STRIP A, $ = 0° , MVD = 16 MICRONS, MASS FLOW = 17 LBM/SEC 

FIGURE 7.30 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR BOEING 737-300 INLET AT a = 15° (^=171 mph, K 0mvd = 0.062) 

(PAGE 6 OF 10). 
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LOCAL IMFIT^OElMCElSrT EFFICIENCY 


TEST RUN ID: 1683615BH, 169361SBH, 1703615BH, 1713615BH, 1723615BH, 1733615BH 



(Positive Surface Distance, S, on Inner Surface) 
fimax = 0.640, So mas = -Hmm, S/ m „ = 44mm ( Experiment ) 

(B2) STRIPS B,H, 0 * 45° , MVD = 16 MICRONS, MASS PLOW = 17 LBM/SEC 

FIGURE 7.30 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR BOEING 737-300 INLET AT a = 15° (^*,=171 mph, K 0m , d = 0.062) 

(PAGE 7 OF 10). 
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LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 1683615CG, 1693615CG, 1703615CG, 1713615CG, 1723615CG, 17336I5CG 



(Positive Surface Distance, S, on Inner Surface) 

0max = 0.610, S 0max = -16 mm, S Imax = 30mm ( Experiment ) 

(B3) STRIPS C,G, 0 = 90° , MVD = 16 MICRONS, MASS FLOW = 17 LBM/SEC 

FIGURE 7.30 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR BOEING 737-300 INLET AT a = 15° (Voc=171 mph, K 0mvd = 0.062) 

(PAGE 8 OF 10). 
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LOCAL IMPINGEMENT EFFICIENCY 


TEST RUN ID: 1683615DF, 1693615DF, 1703615DF, 17I3615DF, 1723615DF, 1733615DF 



(Positive Surface Distance, S, on Inner Surface) 

@m«* = 0.540, S 0m , „ = -55mm, S Jmax = 23mm ( Experiment ) 

(B4) STRIPS D,F, 0 = 135° , MVD = 16 MICRONS, MASS FLOW = 17 LBM/SEC 

FIGURE 7.30 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR BOEING 737-300 INLET AT a = 15° (Voo=l71 mph, K 0mvd = 0.062) 

(PAGE 9 OF 10). 
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LOCAL IMPINGEMENT EEEICIE3STCY 


TEST RUN ID: 1683615E, 1693615E, 1703615E, 1713615E, 1723615E, 1733615E 



(Positive Surface Distance, S, on Inner Surface) 

Pmaz = 0.580, So max = -50 mm, S tmai = 1 6mm ( Experiment ) 

(B5) STRIP E, 0 = 180° , MVD = 16 MICRONS, MASS FLOW = 17 LBM/SEC 

FIGURE 7.30 

AVERAGED LOCAL WATER IMPINGEMENT EFFICIENCY DATA 
FOR BOEING 737-300 INLET AT a = 0° (^=171 mph, K 0m , d = 0.062) 

(PAGE 10 OF 10). 
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8.0 CONCLUSIONS 


Water droplet impingement data were presented for several modern aircraft geometries: 
a Natural Laminar Flow airfoil NLF(1)-0414, two simulated glaze ice shapes, a medium 
speed 30 degrees swept MS-0317 airfoil, a 30 degrees NACA 0012 swept wing tip and a 
Boeing 737-300 engine inlet. 

The work presented in this report is the conclusion of an extensive experimental inves- 
tigation initiated in 1984 to provide experimental water droplet impingement data for the 
validation of computer codes used in icing analyses. The main objectives were to develop 
experimental and data reduction methods for obtaining water droplet impingement char- 
acteristics on aircraft surfaces and to provide experimental impingement data for a range 
of geometries, flow and spray cloud conditions. 

The following conclusions are based on the results obtained during this research pro- 
gram: 

1. The experimental investigation conducted has resulted in an extensive data 
base of water droplet impingement data for several aircraft geometries in- 
cluding engine inlets for which very little data existed prior to this work. 

2. The experimental and data reduction methods are accurate and repeatable. 

The quality of the experimental results depends on careful monitoring of 
all test variables particularly those affecting spray cloud properties. Cloud 
uniformity and MVD size are difficult to control and measure. To improve 
the accuracy of the results each test condition must be repeated several times 
and the results must be averaged to produce the final impingement curves. 

The data reduction method is very efficient compared to colorimetric analysis 
and allows extraction of data without destroying the blotter strips. It is 
important that the laser calibration curve used in the data reduction process 
be checked periodically since the laser characteristics can change over long 
periods. 

3. In water droplet impingement tests accurate LWC measurements are essen- 
tial in obtaining impingement efficiency results. It is important that local 
LWC measurements be made in order to eliminate errors due to variations 
in cloud uniformity. The collector mechanism used in the 1985 tests and 
the improved design used in the 1989 tests proved useful in measuring local 
LWC. 

4. In most cases the uncertainty in the experimental results obtained was less 
than 20%. The maximum uncertainty was in 0 ma x • 

5. Correlation between experimental and computational results was good for 
the majority of cases presented. 
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APPENDIX A - DATA REDUCTION FOR CURVED STRIPS 


Test blotter strips used in the April 1989 droplet impingement tests were of different 
shapes and sizes depending on test geometry. The various blotter strip shapes used are 
shown in Figures A.l, A.2 and A.3. The shape of a blotter strip is important in that 
it affects the data reduction process. The laser reflectometer was designed for straight 
blotter strips only. To scan curved blotter strips each strip had to be divided into straight 
segments and each segment had to be scanned separately. In most cases each curved blotter 
strip was divided into two straight segment. 

Straight blotter strips (Fig. A.l) were the esiest to reduce. These strips were attached 
to the reflectometer drum and the data was extracted by scanning the strip at three 
different horizontal locations as shown in Figure A.l. The three horizontal scans were 
stored in a single file and were averaged to produce a single droplet impingement curve. 
The surface distance along the strip was identical to the surface distance of the geometry 
to which the strip was attached. 

Figure A.2 shows the type of strips used in the testing of the finite swept wing tip. 
These strips were V-shaped to fit the wing section in the streamwise direction. To ex- 
tract the impingement data from the V-shaped strips the strips were scanned along their 
centerline (line A-B-C in Fig. A.2). Section AB of the strip was scanned first, by posi- 
tioning section AB on the reflectometer drum. The results were stored in a data file. Next 
section BC was positioned on the reflectometer drum and was scanned from B to C. The 
impingement results from section BC were stored in a data file. The two data files (one 
for section AB and one for section BC) were combined to produce the final impingement 
curve. To simplify the data reduction process only one horizontal scan was performed for 
each V-shaped blotter strip. 

Figure A.3 shows the type of curved strips used in all impingement tests of the infinite 
swept MS(1)-0317 airfoil. The strip centerline (A-B-C-D-E) is shown in Fig. A.3. The 
centerline is curved so a large number of straight segments were required to preserve the 
actual surface distance. This significantly complicates the data reduction process. A 
simpler and equally accurate method was to approximate the centerline of the strip with 
two straight line segments a-b-c whose total length was equal to the curved centerline A- 
B-C-D-E. The laser scan was performed along the scan line a-b-c using the same procedure 
as the one used for the V-shaped strips discussed above. 
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FIGURE A.l 


STRAIGHT BLOTTER STRIPS USED ON NLF(1)-0414 AIRFOIL, 
ICE SHAPES, BOEING 737-300 ENGINE INLET AND COLLECTOR. 









APPENDIX B - IMPINGEMENT TEST LOG SHEETS 


Details of tunnel and spray system conditions for all impingement test runs performed 
in April of 1989 are given in this appendix. A total of 214 test runs were performed which 
are subdivided as follows: 


• Uniformity tests- Runs 1 to 59 and 181 to 183 

• Boeing 737-300 Inlet - Runs 60 to 71, 93 to 94 and 168 to 180. Runs 60 to 71 
and 93 to 94 were discarded due to problems with the inlet mass flow system 
and the spray system air pressure. 

• Reference Collector - Runs 95 to 115 

• Ice Shapes - Runs 116 to 126 

• NACA 0012 Swept Wing Tip - Runs 127 to 147 

• MS(1)-0317 Swept Airfoil - Runs 148 to 167 

• Laser Calibration - Runs 184 to 192 

• NLF(l)-0414 Airfoil - Runs 215 to 236 

B.l - Test Variables Recorded 

At the top left corner of each log sheet the test date and the dye concentration used 
are given. A dye concentration of 0.00025 was used for the first 59 (uniformity) runs. The 
concentration was then lowered to 0.0002 for the remaining tests. The higher concentration 
allowed easier detection of color variation on the blotter squares during the uniformity 
tests. The lower concentration was selected for the impingement tests to give optimum 
color intensity distribution on the blotter strips. A concentration of 0.0002 was typically 
obtained by combining 4 grams of blue dye powder with 20 liters (i.e., 20,000 grams) of 
distilled water. 

B.1.1 - Run No and Run ID 

Each test performed had a unique number assigned to it. The run number and the run 
ID which was used to identify each model, were placed on each blotter strip used in a test. 


142 



B.1.2 - Tunnel Conditions 

Tunnel conditions, which were measured by the usual IRT instrumentation, were 
recorded for each test point. The tunnel total temperature was obtained from the av- 
erage of eleven readings taken at different tunnel locations. The tunnel true air speed 
along with the freestream static and total pressures were also recorded. The air stream 
humidity level was monitored using the Dew Point temperature reading. A large difference 
between the dew point temperature reading and the air temperature reading indicates a 
low humidity level in the air stream (see table B.l). Two methods were employed to 
increase the humidity in the air stream. In one method moisture was added to the air 
stream by using the IRT spray system. A different approach was to decrease the tunnel 
air temperature by using the IRT refrigeration system. 

B.l. 3 - Spray System Conditions 

Spray system conditions were monitored continuously at various locations in the spray 
system. The three columns under SUPPLY PRESSURE conditions in each log sheet give 
the desired settings for air and dyed water at the spray nozzle. The two columns under 
PLENUM PRESSURE show the spray system pressure settings prior to a test. These 
settings were adjusted manually using pressure regulators located in the tunnel plenum 
chamber. The dyed water pressure was controlled by adjusting the pressure in the tank 
containing the dye/water solution. This pressure was monitored by a pressure gauge 
attached to the dyed water supply tank. The tank pressure was set 6 psig higher that 
the value required at the nozzle to account for the hydrostatic pressure difference between 
the tank level and the nozzle level. The air pressure was set with air flowing through the 
nozzles and was monitored by a pressure gauge attached to the air pressure regulator. 
Both air and water pressure were checked PRIOR TO EACH TEST RUN. To monitor the 
system pressure fluctuations during test runs electronic transducers (Model 204E, Range: 
0 to 250 psi, Setra Systems Inc, Natick, MASS., USA) were attached at various locations 
of the spray system and their readings were monitored from the tunnel control room. 
Variations in the readings of the transducers from their settings prior to each test were 
carefully monitored. The values given in the log sheets were the readings prior to the test. 
Transducers were attached to the supply tank containing the dye and water solution, to 
the dyed water supply line of nozzle 1, 3, 4, 6 and 7 and to the air line of nozzle 11. 

B.1.4 - Mass Flow Conditions 

Mass flow conditions were recorded for all inlet impingement tests. Output voltage 
from the Kurtz probe was recorded for each run. This voltage was converted into pounds 
mass per second using the following equation: 

Mass flow in Ibmj sec = (Volts — 0.11217)/0.13402 

The equation above was obtained by a least squares fit to the mass flow data provided by 
NASA Lewis. 
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B.1.5 - Time Recordings 

For each test point the spray duration in seconds was recorded along with the t ime of 
the day. Spray times of 2.6 and 4.5 seconds were used corresponding to spray clouds with 
MVD’s of 20 and 16 fim respectively. The actual spray time duration is obtained by adding 
0.33 seconds to the spray times shown in the log sheets. This was because the nozzles do 
not stop spraying immediately after the solenoid valves are closed (see Appendix G of Ref. 
10 ). 
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APPENDIX C - NOMENCLATURE 


A f 

Ainlet 

Aoo.m 


C D 

CDinc 

C 

d 

di 
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He — Ne 
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LWCa 

LWCoa 
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MVD 
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n(x) 


n oo 
Pair 
P. »t 

Pt 
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front area of a body, projected parallel to freestream velocity direction 
area of inlet 

area, perpenictxlar to freestream direction, defined by tangent trajecto- 
ries 

chord length 

droplet drag coefficient 

incompressible sphere drag coefficient 

dye concentration 

droplet diameter 

droplet diameter, group i 

infinitesimal impingement area on body surface, Figure 2.2 

infinitesimal freestream area corresponding to dA s , Figure 2.2 

total impingement efficiency in clouds of nonuniform droplet size 

Helium - Neon gas-type laser 

Icing Research Tunnel 

unit vectors Jj to z, y, z axis 

inertia parameter, Equation (2-1) 

modified inertia parameter based on mean volumetric diameter, XK/X a , 
dimensionless 

characteristic dimension of body 

liquid water content at body surface 

liquid water content of droplet cloud at freestream 

Mach number of air flow relative to droplet 

Mean Volumetric Diameter of the droplet cloud 

number of droplet groups available in a discrete droplet distribution 

number of droplets for drops ize diameter x (or d) 

unit normal vector at body surface (pointing outward) 

unit vector parallel to V ^ 

spray system air pressure 

static pressure of air 

total pressure of air 

spray system water pressure 

wall static air pressure 

freestream dynamic pressure 

velocity head experienced by a droplet in the flow field, Equation (D-2) 


R reflectivity 

R too free stream Reynolds number with respect to droplet diameter 

R v Renolds number or air flow relative to droplet 
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VI 

V2 


r oo 

w 

Wp 

Wt 

X 

x,r,0 

x,y,z 

J/oo 


Poo ,m 
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pp 
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distance along surface from point of reference 

distance along surface from point of reference to impingement limit (p), 

Figure 2.2 

True Air Speed 

static air temperature 

total air temperature 

time, dimensionless with L/V, 

time 

droplet velocity, dimensionless with 
droplet velocity 

potential flow velocity dimensionless with Voo 

air velocity at inlet plane 

photodetector voltage output {reflected beam) 

photodetector voltage output (optical feedback) 

freestream airspeed 

freestream air velocity 

engine inlet mass flow 

weighted local impingement rate 

total liquid water content of droplet cloud 

thickness of specimen, Equation (5-3) 

cylindrical coordinates 

Cartesian coordinates 

ordinate of coordinate system perpendicular to freestream velocity di- 
rection at x = —oo, Figure 2.3 
poo corresponding to tangent trajectory. Figure 2.3 
angle of attack 

local impingment efficiency for uniform droplets distribution 

local impingment efficiency for clouds with non-uniform droplet size, 

Equation (2-4, 2-5) 

ratio of specfic heats of air = 1.4 

static pressure difference, Figure 3.7 

Kronecker delta 

central angle of cylinder measured from x-axis 
0 for tangent trajectory 

true range of droplet as projectile injected into still air 

range of droplet as projectile following Stoke’s law 

absolute air viscosity 

density of air 

particle (droplet) density 

density of water 

density ratio 
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